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ABSTEACT 


This report describes programs for the IMAGE-100 system of INPEj 
for extraction of spectral and texticral characteristics ^ followed by 
classification of LANDS AT earth resources multispeatral images. 
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CAPtTULO I 




INTRODUCAO 


1.1 - GENERAL I DADES 


0 campo de processamento digital de imagens tem si do 
objeto d,e muita pesquisa nos ultimos anos. Este fato esta diretamente 
relacionadq ao constante aperfeigoamento dos computadoreS ciigitais, que 
a cada dia se tornam mais rapidos e com maior capacidade de memoria 
principal. Isto tem tornado possTvel a realizagao em tempos, senao cuj2 
tos, pelo menos aceitaveis de di versos tipos de operagoes sobre os da. 
dos de imagem. Entre os objetivos dessas operagoes pode-se citar: 

1) Restauragao ; consiste na recuperagao de uma imagem que tenha 
sido degradada por urn processo qualquer (por exemplo, ruTdo). 

2) Realce : atraves de diversas tecnicas, procura-se melhorar a 
qualidade da imagem ou ressaltar determinados aspectos, a fim 
de facilitar a detegao de certas estruturas ou objetos, elimj^ 
nando informagao irrelevante. 

3) Codificagao : nesta area de atividades, procura-se desenvolver 
tecnicas de representagao de uma imagem, reduzindo o volume 
de dados a serem transmit! dos ou armazenados. 

4) Reconhecimento de Padroes : trata-se do problema da classifies 
gao de uma imagem ou subimagem, em uma as varias classes pre- 
-especificadas 

Talvez pelo fato de ser uma area de estudos bem receji 
te, a maior parte das pesquisas realizadas se limitam a solugao de pr£ 
blemas bem especTficos, envoi vendo certas classes parti cul ares de im£ 
gens, e, so gradualmente, uma general izagao de tecnicas de processameji 
to digital de imagens se form§ (Rosenfeld, 1969). 
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Dentre as iniimeras aplicagoes do estudo de processanieji 
to de ima{|ensj tem-se a destacar aquelas que utllizam imagens iDultie^ 
pectrais, obtidas atraves de sensores remotos a bordo de satelites ou 
avioes. Dentro dos respectivos campos, terfi-se: 

1) Goologia ; Identificagao de tipos de terreno, exploragao de mi_ 
nerais, local izagao de falhas, estudos de erosao. 

2) Agricultura e Agronomia ; identificagao e planejamento de cuj^ 
turas, analise de prejuTzos causados por incendios, geadas e 
pragas, control e de desmatamento , 

3) Ecologia : detegao e monitoramento de poluigao, analise de qu^ 
lidade de agua. 

4) Hidrologia e Hidrografia : monitoramento de cursos de rios, m£ 
peamento da costa. 

5) Geografia : planejamento rural e urbano, planejamento de estr^ 
das, pesquisa de recursos. 

6) Cartografia : construgao e revisao de mapas. 

Pel a importancia e diversidade de aplicagoes, pode-se 
avaliar o impulse dado pelo desenvol vimento de tecnicas de sensoriameji 
to remote a area de processamento digital de imagens. De fato, devido 
a grande quantidade de dados obtidos atraves desses sensores , tornou-se 
necessaria uma utilizagao de meios autoinaticos para processar e anal_i_ 
sar esta informagao. 

Uma imagem multiespectral contem na realidade varies 
imagens (canais) da mesma cena, onde cada canal contem a resposta, em 
termos de reflectancia daquela regiao, a uma faixa estreita do espectro 
eletromagnetico incidente sobre ela. Normalmente, sao utilizados 4 
nais (0.5-0. 6, 0.6-0. 7, 0.7-0. 8, 0.8-l.lym), correspondendo a faixas 
nas regioes de verde, vermelho, infravermelho proximo e intermediario. 
Assim, urn ponto da imagem pode ser interpretado como urn ponto no esp^ 


go espectral 4-dimensional, onde cada coondenada representa o nTvel de 
reflectancia para a cor correspondente a regiao que ele representa. 

A classifica^ao automatica de imagens multiespectrais 
baseia-se no fato de que cada material apresenta uma resposta espectral 
unica. Isto e ilustrado na Figura I.l, onde se tern as respostas tipj_ 
cas da agua, do solo e da vegetagao. Desta forma, a informagao dos qu^ 
tro canais, em conjunto, e utilizada para tentar separar tipos distijn 
tos de materials, de modo que materials diferentes ocupem regioes di£ 
tintas no espago espectral 4-dimensional (assinaturas) . Se as diversas 
classes sao especificadas atraves de assinaturas ou atraves de densida 
des de probabilidade sobre o espago espectral, entao uma classe e atr^ 
buTda a urn ponto, se este pertence a assi nature correspondente aquela 
classe, ou se ela e a classe mais provavel . 

No prime! ro caso, tem-se urn procedimento determinTstj[ 
CO que e o metodo utilizado, originalmente, pelo sistema IMAGE-100. A 
vantagem desse metodo e a possibi 1 idade de treinamento e classificagao 
poderem ser feitos para cada classe (tema), separadamente das demais. 
Embora seja possTvel minimizar a sobreposigao de assinaturas, atraves 
de varios refinamentos na fase de treinamento, se acontecer de urn ponto 
pertencer a mais de uma assinatura, o sistema se omitira da classifies 
gao. No segundo caso (Clessificagao estatistica pelo sistema MAXVER d£ 
desenvolvido no INPE) este problema e superado, obtendo-se precisao 
bastante superior no processo. 

A classificagao, discutida acima, e em geral a mais ut^ 
lizada, principalmente, em imagens multiespectrais. Ela e feita numa 
base de cela de resolugao, uma vez que cada ponto da imagem e classify 
cado isoladamente, nao sendo para isso considerada a vizinhanga do me^ 
mo. 


Pode-se tambem fazer a classificagao numa base de bl^ 
CO de cel as de resolugao contTguas. Urn exemplo seria o caso da identj_ 
ficagao do tipo de culture de urn campo inteiro de plantagao. A chave 



Fig. I.l - Respostas espectrais tTpicas 
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da classiffsagao deste tipo e a definigao de um conjunto de caractem 
ticas significativas, que contenham a i nformagao de interesse, contida 
no bloco. pbviamente este conjunto e muito dependente da aplicagao que 
se tem em vista. Por exemplo, um bom conjunto de caracteristicas, para 
fins de classificagao de microfotografias de celulas do sangue, talvez 
nao fosse viavel para diagnostico automatico de radiografias do torax, 
nem tampouco para reconhecimento do tipo de terreno de certa regiao. 
Porem, uma vez escolhido esse conjunto de caractertsticas, pode-se utj^ 
lizar uma das diversas tecnicas de reconhecimento de padroes, para 
efetuar a classificagao propriamente dita (categorizagao). 

No presente trabalho, so consideraremos a classificagao 
feita numa base de bloco de celas de resolugao contTguas (janela). Por 
tanto, as etapas distintas, em que um sistema classificador desse tj[ 
po se divide, sao: 

Py'ept"ocessamento ; nesta fase, a janela a ser classificada pa^ 
sa por uma ou varias operagoes, de modo a evidenciar determ2 
nadas caracterTsticas que serao usadas para representa-la, ou 
entao, com a finalidade de tornar estas caracterTsticas, inv£ 
ri antes em relagao a posigao, orientagao, etc. 

2) Extragao de caracterTstica s: consiste na realizagao de uma s£ 
rie de medidas sobre a janela preprocessada, com a finalidade 
de se obter um vetor de caracterTsticas ( ou primitivas) un^ 
CO, que represente aquela janela. A finalidade desta etapa e 
a redugao da complexidade do sistema, cujo desempenho vai d£ 
pender de uma boa escolha das caracterTsticas a serem extr^ 
Idas (estas devem caracterizar propriedades comuns a todos os 
membros de cada classe). Como ja foi dito, esta e a parte 
mais dependente daquilo que o sistema se propoe a classificar. 

3) Categorizagao : como resultado da segunda etapa, tem-seumaen^ 
pla de numeros reais, representando a janela em questao. Esta 
enupla pode ser imaginada como um ponto dentro do espago n-dj^ 
mensional. Tem-se, entao, um problema tTpico da area de rec£ 
nhecimento de padroes, onde existe um grande numero de metodos 
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, para o projeto de categorizador (Tou and Gonzales, 1974j Fu 
kunaga, 1972; Meisel, 1972). 

A Figura 1.2 apresenta um sistema pratico para classy 
ficagao automatica de imagens 


IMAGEM 


REAL 


PREPRO- 

CESSADOR 


A 


EXTRATOR 0E 
CARACATERtS 
TICAS 



CATERGO 

y 

r 

RIZADOS 


IMAGEM 

PREPROCESSADA 


VETOR DE 

CARACTERISTICAS 


A 


CATEGORIA 
ATRIBUIDA A 
IMAGEM 


Fig. 1.2 -> Sistema para classificagao automatica de imagens, 


Sendo uma imagem o objeto da classificagao, e natural 
que se utilizem caracterTsticas simi lares aquelas usadas pelo ser hum£ 
no, na interpretagao visual de uma fotografia. Para isto, ele utiliza 
tres propriedades fundamentals: tonal idade, textura e contexto. Estas 
caracterTsticas sempre estao presentes numa imagem, embora eventualmeji 
te uma del as possa estar mais evidente. Da mesma forma, tem-se 3 tipos 
de caracterTsticas que podem ser utilizadas para efeito de classifies 
gao: 

- espectrais: descrevem, em termos estatTsticos , a distribuigao 
dos niveis de cinza dos pontos do bloco considerado. 

- textura: contem informagao sobre a distribuigao especial das v£ 
riagoes de nTveis de cinza, dentro de um canal. 

- contextuais: trazem informagao derivada dos pontos oublocosde 
imagem que estao mais proximos ao bioco que esta sento analj_ 
sado. 

E interessante notar que, embora seja facil para um o^ 
servador humano reconhecer e descrever, em termos empTricos, a textura. 
ela tern si do refrataria a uma definigao precisa e a uma analise exata 
por computadores digitais (Haralick el al., 1973 c). 





No CapTtUiO II, tratar-sa-a da textura e dos tipos 
mais comuns de tecnicas utilizadas para medi-la. 

1.2 - OBJ ET I VOS 

Haralick et al , , (1973 c) sugen'ram uma serie de medj_ 
das de textura, baseadas nas matrlzes de coocorrencia ou de dependeji 
cia espacial de nTveis de cinza. Os resultados encontrados mostraram a 
utilidade dessas caracterTsticas na classificagao de microfotografias, 
de rochas e de tipos de terreno, ero imagem multi espectrais de satelj_ 
tes e em fotografias aereas. 

Este trabalho consiste na implementagao de programas 
no Image-100, visando a extragSo das caracte"5sticas de textura, cita 
das, em conjunto com caracterTsticas espectrais, e subsequente utilize 
po na classificagao de tipos de terreno. A imagem utili.zada foi uma 
do LANDSAT de nP 175152-123730-7, obtitVa na passagem do dia 1/6/75. Urn 
mapa geologico foi considerado como a vs^rdade terrestre, para efeito 
de analise de desempenho das caracterTsitcas . 

Foi desenvolvido tambem i.\m programa para visual izagao, 
no terminal Tektronix, da capacidade de Siwparagao de duas classes quai£ 
quer, utilizando uma determinada caracterVstica. 

0 computador utilizado foi urn POP 11/45, que e o contrjo 
lador de processes do sistema IMAGE-100, coai uma memori a central de 48K, 
A linguagem usada em todos os programas foi FORTRAN- IV. 


CAPlTULO II 


TEXTURA 


2.1 - TEXTURA VISUAL 


0 mundo que nos cerca contem uma grande variedade de 
texturas. Defato, a textura e uma propriedade inerente a todas as supej^ 
fTcies. Pode-se qualifica-la como sendo fina, grossa, irregular, lisa, 
linear e mais urn sem numero de adjetivos. 

Embora todos tenham nogao do significado de textura, e 
muito difTcil definT-la em termos cientTficos precisos. Pickett (1970) 
observe que, para urn padrao visual ser visto como uma textura, ele deve 
ter as seguintes caracterTsticas: 

1) urn grande numero de elementos (variagoes espaciais em intensi_ 
dade ou comprimentos de onda); 

2) os elementos e as regras de espagamento ou arranjo podem ser 
arbitrariamente mani pul ados, desde que uma caracterTstica de 
repetitividade ainda continue; 

3) desde que haja suficiente detalhe em urn pequeno angulo de vis^ 
da, uma textura caracterTstica desponta, mesmo quando os elemeji 
tos basicos ou os espagamentos sao aleatoriamente distribuTdos. 

Ve-se que os atributos basicos de uma textura visual 
sao muitas e vepetitivas vaviagdes 

Nao e ainda bem conhecida a maneira pel a qual o ser 
mano deteta e analisa informagao textural, no processo de percepgao vj[ 
sual de uma cena. Julesz (1965 e 1975) considera que existem dois n^ 
veis distintos de percepgao visual de textura. Urn espontaneo ou im 
pressionista, caracterizado por uma percepgao pura e outro deliberado, 
que exige a utilizagao de processo cognitivos para ser completa a per 
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cepgao. Trabalhando com imagens artificials obtidas por computador, 
ele investiga a importancia da textura na percepgao visual Humana e an£ 
lisa as diversas limitagoes dessa percepgao. 

Tambem na aplicagao especTfica em imagens, nae se encojv 
tra urn metodo perfeito de analise de textura. Sabe-se, entretanto, que 
el a desempenha papel importante, na detegao e no reconhecimento de o^ 
jetos (Pickett, 1970), e na discriminagao de tipos diferentes (Har^ 
lid et al . , 1973 c e Haralick e Bosley, 1973 a). Hawkins (1970) mostra 
as dificuldades existentes no trabalho de extragao de informagao textu^ 
ral e conclui que a classificiacjo de texturas e uma das mais dificeis 
tarefas, no campo de processamehto de imagens. 

A textura pode ser analisada em dois niveis: o estate 
tico e 0 estrutural . Do ponto de vista estatTstico, a textura e define 
da atraves de urn conjunto de parametros estatrsticos, obtidos deumgra^ 
de numero de medidas locais, feitas na imagem. Do ponto de vista estrjj 
tural, tenta-se local izar elementos que ocorrem repetidamente e desc^ 
brir as regras de organizagao dos mesmos, dentro de determinadas suba^ 
reas da imagem. Quando se trabalha com imagens naturals, geralmente se 
utiliza 0 metodo estatistico, pois o estrutural e mais complicado. 

A seguir serao mostradas algumas das tecnicas mais c^ 
muns, usadas para fins de classificagao de imagens, no sentido consid^ 
rado neste trabalho, que e o de discriminagao de blocos de celas de r^ 
solugao, em categories preespecificadas. 

2.2 - ESPECTRO E POTENCIA DE FOURIER 


A transformada de Fourier, de uma imagem f(x,y), e de 


f ini da como: 


•f-OO^-OO 


F(.,v) . e ■ (''X + vy) 
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0 valor da potencia num ponto (u,v) e: 

|F(u,v)|^ = F(u,v) . F* (u,v) 

As caracteristicas de textura, baseadas em medidas do 
espectro de potencia de Fourier, utilizam dois fatos: 

1) A distribuigao radial ( no espago tranformado) dos valores de 
potencia esta relacionada a rugosidade da textura. Enquanto 
uma textura grossa tern valores altos de potencia perto da ori_ 
gem, numa textura fina, os valores de potencia sao mais esp£ 
lhados. 

2) A distribuigao angular dos valores de potencia esta relaciona 
da a direcionalidade (existencia de bordas e linhas numa ce£ 
ta diregao) da textura. Uma textura orientada numa diregao 0 
tera valores altos de potencia, concentrados perto da perpeji 
dicular a esta diregao. 

Este fatos sugerem que sejam usadas caracteristicas de 
textura, que sejam relacionadas as medias da potencia, ao longo de 
aneis centrados na origem (para analise de rugosidade) e ao longo de 
cunhas (para analise de direcionalidade) 

A transformada discreta de Fourier, deumaimagem dig_i_ 
tal de N X N pontos, e: 

F(u,v) = — V f(m,n) e 

m=0 n=0 


0 ^ u,v $ N-1 










Tem-se, entao, dois tipos padroes de medida de textura, 
baseados no espectro de potencia: 


, r = — ' |F(u,v)P, r e r arbitrarios. 

1^2 ^ jj2 y2 ^ yi2 1 2 

1 2 

- ■ |F{u,v)|^, 0 e0 arbitrarios. 

*0,0 =0 tg" (v/u) < 0 ^ ^ 

Pode-se tambim fazer uma combi nagao desses tipos, oque 
corresponderia a calcular o somatorio dos valores de potencia, na inte£ 
segao de aneis e cunhas. 

Pelo fa to de, em geral , a imagem nao ser periodica, seji 
do seus lados opostos descontTnuos, aparece uma direcionalidade adicio 
nal ra horizontal e na vertical, que nao deve ser levada em conta. 

2.3 - ESTATTSTICAS DE DIFERENQ,-.* DE NTVEIS DE CINZA 

Considere-se uma imagem digital f(m,n) quantizada em N 
nTveis de cinza. Para urn deslocamento arbitrario d= (Am, An), onde Am 
e An sao inteiros, e computado urn vetor N-dimensional p^, da seguinte 
forma: sua i-esima componente p^ (i) e o numero de pares de pontos da 
imagem separados de d, tendo uma diferenga de nTveis de cinza igual a 
i. Se este vetor for normalizado, p^- (i) fornece uma estimativa da pr^ 
babilidadede cue fjj(m,n) = |f{m,n) - f(m + Am, n + An) | assuma o valor i. 

Observe-se que, para uma textura grossa e d pequeno, 
comparado ao tamanho do elemento de textura, e esperado que os pares 
de pontos separados de d tenham, normalmente, nTveis de cinza proximos, 
de modo que p^ (i) assume valores altos para i's baixos. Da mesma for 
ma, se a textura for fina e d for comparavel ao tamanho do elemento, 
entao, espera-se uma melhor distribuigao dos valores das componentes 
de Pjj. 
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Medindo-se o grau de espalhamento dos valores em p^ 
ra varios d's pode-se, entao, ter uma medida da rugosidade da textura. 
tem-se por exemplo, as caracterTsticas: 

N-1 

- Contraste: CON = E p.(i) 

i=1 ° 

, N-1 

- Media : MED = Z i . Pj(i) 

N i=1 ° 

Alem dessas caracterTsticas» pode-se obter informa?ao 
adicional sobre a direcional idade da textura, quando se comparam os dj_ 
versos graus de espalhamento de p^ a medida que a diregao de d e vari£ 
da e sua magnitude e mantida constante. 

Outro conjunto de caracterTstica deste tipo pode ser 
obtido atraves de vetores p^ extraTdos, nao a partir de pares de pontos 
isolados, mas fazendo-se uma media de nTveis de cinza sobre urn grupo de 
pontos vizinhos. 



2.4 - ESTATTSTICAS DE COMPRIMENTOS DE CADE IAS 


I 

^ I 


Dividem-se os N nTveis de cinza em M faixas distintas. 
Uma matriz e, entao, obtida de modo que seu elemento p(i,j) seja o nQ 
mero de cadeias de comprimento j, em uma diregao 0, consistindo de poji 
tos, cujos nTveis de cinza estejam dentro da i-esima faixa. 


E bem intituivo o significado dessas cadeias: se uma 
textura e grossa, espera-se encontrar cadeias longas, ao passo que se 
a textura for fina, as cadeias serao mais curtas. 

Util izando-se parametros convenientes dessa matriz (por 
ex(?mp1o, urn que dependa fortemente do numero de cadeias longas), e pos 
sivel obter medidas do tipo de textura da imagem. A direcional idade p£ 
de ser medida atraves do levantamento de varias dessas matrizes, uma 
para cada diregao 0 das cadeias. 






('.•a;/- 





CAPlTULO III 


ESTATlSTICAS DE OCORRENCIAS 


No segundo capTtulo foram citados dois tipos de medidas 
de textura, baseadas em estatisticas de propriedades locals. Num-, estas 
propriedades eram diferengas de nTveis de cinza e no outro, comprimeri 
to de cadelas. 


Neste capTtulo, e mostr^ido com mais detalhes urn tercej^ 
ro tipo de estatisticas proposto por Haralick (1973 c), baseado nas 
trizes de dependencia espacial de nTveis de cinza, ou matrizes de c£ 
ocorrencia. Assume-se que toda a informagao de textura, existente, pode 
ser convenientemente representada (ou codificada) por urn conjunto de 
matrizes computadas, para varias diregoes e distancias entre pares de 
cel as de resol ugao da imagem. 

3.1 - MATRIZES DE COOCORRENCIA 


Suponha-se que a imagem digital em analise, com N ni] 

veis de cinza, seja retangular e const! tuTda de N cel as de resol ugao 

A 

na diregao horizontal e N^ cel as de resol ugao na diregao vertical. Uma 
tal imagem sera denominada janela. Aqui , nao se considera pontos, mas 
cel as de resol ugao, porque pode-se tambem trabalhar com vizinhangas de 
pontos elementares, e portavVto, com nTveis de cinza medios nestas viz_i_ 
nhagas. Sejam: 

L = {1,2, ..., N } ... domTnio espacial horizontal 

X X 

L = {1,2, ..., N } ... domTnio espacial vertical 

6 = 1,2, ..., N ... conjunto dos nTveis de cinza possiveis. 

L X L ... conjunto de todas as cel as de resol ugao da janela. 

y ^ 
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Nesta notagao, a imagem seHa uma ftmgao I,queassocia 

a cada cel a de resol ugao, algum nml de cinza, isto e, I:L x L -*■ G. 

X y 

Cada matriz e construida de tnodo que, para utna diregao 
0 e uma distancia especificada d, seu elemento P (i.j) seja a freque_n 
cia de ocorrencia na janela de pares de resolugao distanciados de d, 
numa posigad a 0 °, uma com nTvel de cinza i e a outra com nTvel de cin 
za j . 


Sao usadas diregoes 0 iguais a 0°, 45°, 90° e 135°. De 


finindo formal mente cada matriz, tem-se: 


P(U.d,0°) = # {((k,l,),(m.n)) e (L^ x ^ x (L^ x L^) 
I(K,A) = i, I(m,n) = j} 


k-m= 0 , [l-n|=d. 


P(i,j,d,90°) = #= {((k,l,),(m,n))e(Ly x 4 ) x (L^ x L^) 
I(k,l) = i , I(m,n) = j} 


k-m|=d,l - n= 0 , 


P(iJ,d,45°) = # {((k,l),(m.n) e (L^ x L^) x (L^ x L^) 


(k-m=d.l-n=-d) 


ou (k-m = -d,l-n=d), I (k,l) = i, I (m.n) = j} 


P{i,j,d,135°) = # C((k,l,),(m,n)) e (Ly X y x (Ly x l^) (k-m=d,l-n=-d) 
ou (k-m=-d, l-n=-d), I (k,l)=i, I(m,n) = j} 


onde # signiflca 0 numero de elementos do conjunto. 

Na realidade, como se pode ver das defi ni goes das matri^ 
zes, 0 conceito de distancia, utilizado, nao foi euclideano: 


dist ((k,l ) ,(m,n)) = max {|k-m|,|k-nl> 
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Pode-se tambem observar que essas matrizes sao simetH 
cas, istq e, P(i, 3 ,d, 8 ) = P(j,i,d, 0 ). 

A Figura III.l mostra uma imagem simples 4 x 4, com 3 
nTveis, e as matrizes correspondentes no caso da distancia ser unitaria. 


112 3 

4 3 1 


2 1 1 

1 2 3 3 pQ = 

3 4 3 

*^45 " 

1 4 3 

2 2 2 3 

1 3 2 


1 3 2 

3 1 1 2 





2 4 o] 


0 4 1 

^90 "" 

4 2 3 

^135 

4 4 1 


0 3 4 


1 1 2 


Fig. III.l - Imagem simples e matrizes de coocorrencia 
correspondentes. 

Na determinagao, por exemplo, do elemento (3,2) da m^ 
triz horizontal (8 = 0 °), conta-se o numero de pares de celas adjanceri 
tes, horizontalmente, tal que a primeira tenha nTvel 3 e a segunda nj_ 
vel 2 . 


Devi do a definigao das matrizes , nota-se que cada par 
de celas contTguas e contado duas vezes. A primeira no levantamento da 
entrada (i,j) da matriz, e a segunda no da entrada (j,i), onde i e j 
sao os nTveis de cinza das duas celas. 


Pode-se fazer uma normal izagao de cada uma dessas m^ 
trizes, de modo que a soma de todos os seus elementos unitaria. Para 
isto e necessario o calculo do numero R de celas de resol ugao vizinhas 
(de acordo com 6 e d). No caso de distancia 1, tem-se; 


8=0: para cada linha existem 2(N -1) pares de celas vizinhas, 

X 

Como existem linhas, ha urn total de _i) 
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0 = 90°: analogamente, Ry = 2N^ (N^ -1) 

0 = 45°: para cada 1 i nha , a excegao da primeira, ex1tem2 (N -1 )pa 

A 

res de cel as vizinhas. Entao = 2(N^“l)(Ny-l) 

0 = 135: analogamente, R^ = 2(N^ - ^)(^x ” 

Para uma distancia d qualquer, obviamente roenor que N 

A 

e N , tem-se: 

R^ = 2Ny(N^ - d) Rj 3 = 2(N^ - d) (N^ - d) 

Ry = 2N^(Ny - d) Rg = 2(N^ - d) (N^ - d) 

Apos a normal izagao, pode-se encarar urn certo elemento 
P(i,j) de cada matriz, comO sendo uma estimativa da probalidade de que 
urn dado par de cel as de resol ugao, satisfazendo a relagao espacial para 
a qual a matriz foi computada, tenha nTveis de cinza i e j. Embora a 
normal izagao nao traga informagao adicional sobre a textura da imagem, 
uma vez que os elementos p (i,j) sao iguais aos P(i,j) da matriz nao 
normal izada, a menos que uma constante multi pi icativa adequada, el a se 
apresenta muito util em dois aspectos: 

1) aspecto computacional : apos a extragao das matrizes, e feito 
urn intenso processamento, visando a obtengao de caracteristi_ 
cas de textura. Algumas operagoes envoi vi das podem levarauma 
ultrapassagem da capacidade de representagao de numeros reais 
no computador. Isto ocorre, principalmente, quando a janela e 
constituTda de muitas cel as de resol ugao: 

2) de certa forma a normal izagao torna as caracterTsticas indepeji 
dentes das dimensoes da janela que representam, dependendo ap^ 
nas do tipo de textura contida. Ou seja, duas janelas com tex 
turas identicas, uma estando contida na outra, forneceriam as 
mesmas caracterTsticas. Isto, entretanto, ocorre dentro de cer 
tos limites, pois se a janela menor for muito pequena, ela po 


de nao center informagao de textura suficiente, de mode que 
sen' a mal caracterizada. 

A esta altura e interessante citar dois aspectos rela 
cionadol a (jl^tengao das matrizes de coocorrencia. 0 primeiro e o fato 
de que o numero de operagoes necessarias, para processar uma janela, e 
diretamente proporcional ao numero n de cel as de resol ugao existentes. 
Em comparagao, para se obter transformadas de Fourier ou Hadamard, pa 
ra fins de medida de textura, este numero seria proporcional a n.logn. 
0 segundo aspecto e que so ha necessidade de trazer, a memoria do com 
putador, duas linhas da janela ao mesmo tempo, nao havendo, porta nto , 
forte exigencia de armazenamento neste tipo de processamento. 

3.2 - A EQUALIZACAO DO HISTOGRAMA 


No processo de digitalizagao da imagem, cada valor do 
grau de cinza e quantizado sobre uma faixa finita de niveis possTveis. 
Em consequencia, pode ser computado urn histograma de niveis de cinza 
que, a menos de uma constante, seria uma estimativa da fungao densid^ 
de de probabilidade da variavel nTvel de cinza para aquela imagem. 

A equal izagao do histograma consiste num processo apos 
0 qual , todos os niveis de cinza, da imagem transformada, apresentam a 
mesma frequencia de ocorrencia. 

Embora a textura e a tonal idade sejam duas propriedai 
de que se interrelacionam, e natural e desejavel que duas imagens de 
uma mesma cena, uma com nTveis de cinza iguais aos da outra,amenos de 
uma constante aditiva ( o que equivale a urn histograma deslocado), foj2 
negam caracteristicas de textura iguais. A equalizagao do histograma 
garante que duas imagens, uma sendo transformagao monotonica da outra, 
fornecem os mesmos resultados (Haralick et al . , 1973 c). 

Na pratica, esta transformagao monotonica pode ser d£ 
vida a diferengas de iluminagao, revel agao, filmes, lentes ou digital^ 
zador. 
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Alem disso, a equal izagao permite que se fixe o numero 
de nivets de cinza com que ss esta interessado em trabalhar. Isto e im 
portante, devido a necessidade de armazenamento das matrizes. Existin 
do N nTvels de cinza, e necessario o uso de palavras reais de mem£ 
ria para cada matriz (ou N(N+1)/2, no caso de se utilizer a propried^ 
de de simetria). Acontece que a maioria das imagens apresentam urn hi^ 
tograma concentrado numa faixa estreita de nTveis de cinza. Aqui, nao 
se tratando de acentuar ou melhorar visualmente o constrate, pode-se 
simplesmente trabalhar com urn numero de nTveis aceitavel pelo computa 
dor e ipatTvei com esta faixa estreita do histograma. 

No caso especTfico deste trabalho, a equalizagao adm^ 
te que a imagem tenha no maximo 256 nTveis (resolugao maxima do si sterna 
I-lOO), transformando^se em 32 nTveis equal izados. 0 mecanismo e simj_ 
lar ao uso na transformagao de distribuigoes em EstatTstica. 

Seja 0 histograma da imagem no canal, em que a textura 
vai ser utilizada para efeito de classificagao A(i), i =0,255. A parti r 
dele, obtem-se a fungao de distribuigao acumulada dos nTveis de cinza, 
conforme a Figura 1 1 1. 2. 


i 
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Entao, e geTada uma tabela, relacionando os niveis, na 
imageni original, aos correspondentes apos a equalizagao, de modo que 


32 NT 


s (-^) 

32 


( 1 ) 


entao, ao nTvel i sera assinalado o nTvel 

Pelo fato da distribuigao ser discrete, aparecem pro 
blemas nas transagoes entre niveis da imagem transformada, quando e]i 
tao um certo niveT, na imagem original, teria que ser designado para 
dois niveis diferentes. Isto seria resol vido com o uso, na fase de e 
qualizagao propriamente dita, de uma sub-rotina, de numeros aleatorios 
de mode a se "rachar" um certo nivel em dois outros, fornecendo uma e 
qualizagao mais refinada. Como isto, forgosamente, demandaria um tempo 
cdnsi deravel , resolveu-se fazer, simplesmente, uma escolha pela melhor 
colocagao daquele nivel original em um dos dois niveis transformados , 
nao sendo portando, a equagao (1) seguida a risca. 

Seria interessante citar que, pelo menos na imagem uti_ 
lizada, normalmente uma janela de 62 x 62 pontos nao apresenta mais que 
40 niveis de cinza, com frequencia de ocorrencia nao nula, nao sendo, 
portanto, necessaria a utilizagao de um numero maior de niveis na jan<e 
la equal izada. 


3.3 - CARACTERISTICAS DE TEXTURA 


A seguir, mostra-se como sao obtidas as medidas de te2< 
tura, a partir de cada uma das matrizes de coocorrenci a (Haralick, 
1973c). Como estas matrizes sao simetricas, na definigao de algumas 
destas medidas, este fato ja foi levado em conta. 
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Sejam: 

N... ntimero de nTveis de cinza apos a equal izagao da janela 
p(i,j) ... elemento (i,j) da matriz de coocormncia ja normalizada. 


I 

i 

<:tj 


N 

ALFA (i) = S p (i,j) ... soma dos elementos d?. i-esiift liriha (ou coluna) 

j < 


BETA (k) = 


GAMA (k) = 


N N 

2 2 P(i»d). k 

i=1 j=1 

i+j=k 


N N 

2 E p(i,j), K 

i=l j=l 


i-j =k 


2,3 .... 2N soma dos elementos da k-es^ 
ma diagonal secundaria 


0,1 .,., N-1 duas vezes a soma dos el£ 
mentos da k-esima diagonal. 


T 

.1 

I 


1 


N N N 

u = E E i p(i , j) = E i ALFA (i ) 
i=j j=l i=j 



» 1 

■ N N 

— 

1/2 

" N 

V = 

i 

1 

E E i^ P(i.j) 

-u" 

= 

E i^ ALFA (i) - 


! 

L d=i J 





1/2 


Tem-se, entao, as caracterTsticas; 
1 ) Momento angular de 2a. ordem : (ASM) 


N N 

FI = E E {p (i,j)}=^ 

i=l j=l 


2) Contras te: (CONT) 


N-1 

F2 = E n^ GAMA (n) 
n=0 




3) Corre1acao ;(CQk) 


i = 1 j = 1 i j p ( 


iJl " U^ 


4) Variancia ; (VAR) 


N N 

F4 = E E (i-u)2 p(i,j) 

i=l j=T 


5) Momento inverso da diferenga: (XIDM) 


N N .X 

F5 = E E — 

1=1 j=i 1 + (i-j)2 


6) Me dia da soma ; (SUMAVE) 


F6 = E i BETA (i) 
i=2 


7) Variancia da soma: (SUMVAR) 


F7 = E (i-F6)2 beta (t) 
i=2 


8) Entropia da soma: (SUMENT) 


F8 = - E BETA (i) In (BETA (i)) 






\ "" 
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9) Entropia ; (ENT) 


N N 

F9 = - S L p(i,j) ln(p(i,j)) 

i=1 j=l 


10) Varlancia da Diferenga ; 


N-1 N-1 

no - - E (i-DIF)2 GAMA(i), ondeOIF= E iGAMA(i) 
i=0 i=0 


11) Entropia da Diferenga ; (DIFERENT) 


N-1 

FVl = - E GAMA (i) In (GAMA(i)) 
1=0 


12) Medida adicional de correlagao: (CORME) 


FI 2 = F9 - HXY onde: 

HX 

N 

HX = - E ALFA (1) In (ALFA (i)) 

1=1 

N N 

HXY = - E E \r\{l^^{i) l\L^ 

1=1 j=l 

Alein destas, foram sugeridas mais duas outras medidas 

de correlagao: 


FI3 = (1 - exp {-1(HXY2 - F9)}, onde 









^0 




' i 

^ i 


I 

i 

1 

j 





HXY2 = 


N 

- t 
U\ 


N 

S 

j-1 


ALFA (i) ALFA(i) ln{ALFA(1) ALFA(j)} 



F14 p (29 maior autoval<^r de Q) , 


onde: 





ALFA (i) ALFA (j) 


Iriiclalmente todas as medldas foram implementadas. 0 F13 
geralmente causavA problema de "overflow", mesmo com a normal izagao 
das matrizes e com^tilizacao de logarltmo decimal no calculo de F9 e 
HXY2. 0 FI 4, alOT de exigir uma memoria extra, pels a matriz Q deveria 
ser caiculada com precisao extendi da , tomava muito tempo da maqu1na> 
Por estes motivos nao foram raais utilizadas estas duas medldas, na im 
plementagao finaT do programa extrator de caracterTsticas. Em compens£ 
gao devido a relativa folga de memoria, foi possTvel uma maior intera 
gao entre o usuario e o brograma, alem de grande economia de tempo. 


‘ Embora seja difTcil explicar o que representa exatamente 
cada uma destas medidas, e intuitivo, pelo menos para algumas, o seu 
significado. 


Ao se analtsar o ASM, ve-se que ele e pequeno quando t£ 
dos os elementos da matriz sao aproximadamente iguais, aumentando a m£ 
dida que alguns elementos fi cam mai ores que os outros. Isto pode ser 
encarado como uma medida de homegeneidade, pois nu^lia imagem homogenea 
exitem poucas transigoes de nTveis de cinza dominantes. 

0 CONT mede 0 momento de inercia da matriz em relagao a 
sua diagonal principal. Sendo a medida mais natural do espalhamento dos 
valores da matri z, Note'^se que, se a textura for grossa e a di stand a 
d pequena, emf relagao ao tamanho do elemento de textura, os pares de 
pontos separados de d devem normal men te ter valores proximos de nTveis 

^ ^ V ■■ ■ ■■ " 

■?) 

f 


Cr 
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de cinza. Isto implica em que os elementos, proximos a diagonal princi^ 
pal, deverao assumir valores elevados. Por outro lado, se a textura for 
fina e a distancia d comparavel ao tamanho do elemento da textura, os 
pares de pontos separados de d terao nTveis de cinza mais separados, 
ou seja, havera um espalhamento uniforme dos valores da matriz. 0 CONT 
medira, entao, a quantidade de variagoes locals de nTveis de cinza, pr£ 
sentes na imagem, ou seja, o contraste. 




i / ^ 
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A correlagao (COR) fornece uma medida de quao similares 
sao as linhas ou colunas da matriz de coocorrencia. Se os valores sao 
bem distribuTdos, tem-se alta correlagao, enquanto que, no caso contrS 
rio,. tem-se baixa correlagao (por exemplo, quando os valores longe da 
diagonal sao muito pequenqs). Em termos de textura, a correlagao for 
nece uma medida das dependencias lineares dos nTveis de cinza na ime 
gem ( Haralick, 1973 c). 

Cada uma das medidas citadas e extraTda de uma matriz, 
sendo, portanto, fungao da diregao e da distancia para as quaiselafoi 
obtida. E claro que e bastante informative o fa to de uma certa textura 
ser direcional, apresentando estruturas em uma certa orientagao. Contjj 
do, nao se esta interessado na diregao parti culardessas estruturas. E 
interessante, portanto , que sejam utilizadas caracterTsticas de textu 
ra que sejam invariantes a orientagao da imagem, mas ao mesmo tempo tr^ 
gam informagao de direcional idade. Obviamente, nao seria o caso de se 
utilizer diretamente as caracterTsticas obtidas das quatro matrizes. 

Foram utilizadas fungoes dessas caracterTsticas, da s^ 
guinte forma: supondo que uma determinada janela tenha caracteristicas 
a, b, c e d, correspondendo as quatro diregoes, obteve-se a media, a 
faixa de variagoes e o desvio padrao deste conjunto, os quais sao inv^ 
ri antes a uma rotagao da janela. 

Tem-se, entao, um conjunto de 12 x 3 caracterTsticas de 
textura que agora poderao ser utilizadas como dados de entrada para o 
classificador. E certo que existe uma forte dependencia entre varias 
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dessas 36 medidas, e a utiliza^ao de todas ao mesmo tempo trarla multa 
informagao superflua ao classificador, sendo, portanto, necessaria uma 
selegao de caractensticas. Isto se constituti num problema basico da 
area de reconhedmento de padroes. 

Naoconstado presente trabalho urn programa de selegao 
de caracterTsticas. Entretanto, foi implementado urn programa (HGRAM) 
que permite uma comparagao das distribuigoes de duas classes arbitra^ 
rias, em relagao a uma determinada caracteristica, atraves de uma di£ 
cretizagao dos valores assumidos por esta caracteristica e posterior 
levantamento dos histogramas relativos as classes escolhidas. 


CAPITULO IV 


ESTUDO PRINCIPAL 


4.1 - DADOS UTILIZADOS 


Para este estudo foi utilizada uma imagem do LANDSAT denO 
mero 175152-123730-7, do nordeste de Mato Grosso, obtida na passagem do 
dia 1 de junho de 1975. Um mapa geologico do Departamento Nacional da Pro 
dugao Mineral, do Minist'erio das Minas e Energia, foi usado para separar 
tres tipos principals de terreno existentes: 

Tipo A: Pre-cambriano 
Tipo B: Coberturas detrTticas 
Tipo C: Pleistocene 

A Figura IV. 1 mostra a banda n9 7 da imagem (infravermelho 
intermediario), que foi escolhida para analise de textura. Foram selecio 
nadas 182 janelas de 62 x 62 pontos, sendo 57 da categoria A, 62 da B e 
62 da C. 

4.2 - PROCEDIMENTO 


Usou-se a mesma resol ugao obtida com o LANDSAT, ou seja, 
cada ponto da imagem no sistema I-lOO, correspondente a exatamente um pon^ 
to da imagem do LANDSAT. Como esta e composta de 3240 x 2340 pontos, e a 
capacidade maxima da memoria de imagem do sistema I-lOO e de 512 x 512 
pontos, foram necessarias varias etapas para que toda a imagem fosse pr£ 
cessada. 


Considerando-se a resolugao da imagem do LANDSAT, onde c^ 
da ponto representa uma area de cerca de 70 m x 60m, foi- usada a distan 
cia unitaria para obtengao das matrizes de coocorrencia. 
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Fig. IV. 1 - Imagem LANDSAT do nordeste do Mato (irosso. 
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dualmente da o melhor resultado e, em seguida, uma terceira, que junto 
com aquele par, da o melhor resultado e assim por diante. 

No caso deste trabalho nao se trata do problema da procjj 
ra de urn conjunto otimo de caracterTsticas, mas de uma analise da impor 
tancia das caracterTsticas de textura em geral . Portanto, para uma cla^ 
sificagao com pares de caracterTsticas, foram escolhidas as medias de to 
dos os tipos de medidas de textura (12 no total), alem da faixa de vari^ 
gao e do desvio padrao do XIDM, que foram exatamente as caracterTsticas 
de textura de melhor desempenho individual. Das caracterTsti cas espe£ 
trais, foram desprezadas as duas que deram piores resuTtados, ou seja 
4-V e 5-V. 


Fica-se, entao, com urn total de 20 caracterTsticas para se 
rem utilizadas em classificagoes, atraves de pares de caracterTsticas . Os 
resultados obtidos, quando o classificador usado e o MTn-Max, encontram- 
se na label a IV. 2. 

0 melhor resultado, encontrado com a utilizagao de pares 
de caracterTsticas de textura, foi de 85 acertos em 91 possTveis, ou s£ 
ja, uma prrcentagem de classificagao correta de 93,4^. Is to foi obtido 
com o uso das caracterTsticas XIDM-D e CORME-M. 

Ja 0 melhor resultado, com caracterTsticas espectrais.foi 
de 87 acertos (95,6%), obtido com o uso das caracterTsticas 7-M e 7-V. 

0 melhor resultado geral foi 90 acertos, correspondendo a 
98,9% de classificagao correta, e obtido com urn par misto de caracterT^ 
ticas, ou seja, uma textural (SUMENT-M) e outra espectral (7-M). 

De forma global, os melhores resultados foram os obtidos 
com pares mistos de caracterTsticas. Isto pode ser melhor conipreendi do com 
a ajuda das Figuras IV. 2 - IV. 6, onde estao colocados alguns histogramas 
obtidos, usando-se o conjunto inteiro das 182 janelas. Nestas Figuras , as 
elapses 1 , 2 e 3 sao usadas para designer, respecti vamente»os tipos de terre 
no A, B e C. 






Fig. IV. 3 - CaracterTstica; XIDM-D 










Seja 0 histograma da Imagem no canal, em que a textura 
vai ser utilizada para efeito de classiftcagao A(i), i=0,255. A partir 
dele, obtem-se a fungao de dlstrlbuigao acumulada dos nTveis de cinza, 
conforme a Figura II I. 2. 
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quallzagao mals reflnada. Como Isto, forgosamente, demandaria urn tempo 
consideravel , resolveu-se fazer, slmplesmente, uma escolha pela melhor 
colocagao daquele nTvel original em urn dos dois nTveis transformados, 
nao sendo portando, a equagao (1) seguida a risca. 

Seri a interessante citar que, pelo menos na imagem uti 
Tizada, normalmente uma janela de 62 x 62 pontos nao apresentamais que 
40 nTveis de cinza, com frequincia de ocorrencia nao nula, nao sendo, 
portanto, necessaria a utilizagao de urn numero maior de nTveis na jane 
la equal izada. 

3.3 - CARACTERISTICAS DE TEXTURA 


A seguir, mostra-se como sao obtidas as medidas de tex 
tura, a partir de cada uma das matrizes de coocorrencia (Haralick, 
1973 c). Como estas matrizes sao simetricas, na definigao de algumas 
destas medidas, este fato ja foi levado em conta. 
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A Figura IV. 2 mostra os histogramas para o ASM-M. Pode-se 
ver que existe boa separagao entre as classes 1 e 2, e deduzir-se que a 
separagao e melhor ainda para as classes 1 e 3, Como o terreno do tipo 2 
e mais homogeneo que o do tipo 1, era esperado que os valores assumidos 
pelo ASM fossem mai ores para o tipo 2, o que de fato aconteceu. 

A Figura IV. 3 mostra os histogramas para o XIDM-D, onde se 
pode ver que, embora haja certa intersegao entre as classes 2 e 1, e 2 e 
3, para as classes 1 e 3 ela ja e menor. Na realidade, a existencia de i£ 
tersegao nao implica na inutilidade da caracterTstica, pois o simples fa 
to das medias dos valores desta caracterTstica seremdiferentes, para cla^ 
ses dife rentes, ja mostra que ela e informativa. 

Da mesma forma, podem ser interpretadas as Figuras IV. 4 p^ 
ra 0 ENT-M, e IV. 5 para o CORME-M. 

A Figura IV. 6 mostra as distribuigoes das classes, quando 
a caracterTstica usada e a media espectral da janela na banda 7 (7-M), 
que foi utilizada para extragao de informagao sobre a textura. Pode-se 
ver que a classe 3 e bem separada da classe 2, e deduzir que a separagao 
e um pouco menor com relagao a classe 1. Isto e devido ao fato de que a 
maior parte do terreno tipo C e bem escura, diferenciando-se dos demais 
tipos, 0 que pode ser visualizado na Figura IV. 1. 

Desta forma, ja eram esperados os resultados obtidos, uma 
vez que em geral as caracterTsticas de textura, usadas individualmente, 
separam bem as classes 1-2 e 1-3, enquanto que as espectrais (emespecial 
7-M) separaram bem as classes 2-3. 

As medidas espectrais do tipo variancia, via de regra, ti_ 
veram desempenho mais fraco que as do tipo media, quando usadas em pares 
mistos de caracterTsticas. Porem, notou-se que os pares de caracterTsti_ 
cas espectrais, com melhores resultados, foram aqueles que usaram uma do 
tipo media e a outra do tipo variancia. 
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Em seguida, foi realizado urn estudo de classificagao, uti_ 
lizando-se a dimensao do espago de caracterTsticas igual a 3. .Utilizando- 
se as 20 caracterTsticas anteriores, podiam ser feitas 1140 combinagoes 
diferentes para classificagao. No entanto, foram consideradas apenas 12 
caracterTsticas, sendo 9 texturais e 3 espectrais. As espectrais foram 
7-M, 7-V e 6-M. 0 criterio usado para selegao das 9 caracterTsticas de 
textura t:oi a soma total de acertos dos pares de caracterTsticas, dos 
quais cada uma dessas 9 caracterTsticas fez parte, o que pode ser obtido 
da Tabela IV. 2. Foram, entao, escolhidas as caracterTsticas ASM-M,VAR-M, 
XIDM-M, XIDM-D, XIDM-F, SUMVAR-M, SUMENT-M, ENT-M e CORME-M. As tabelas 
IV. 3 - IV. 11 mostram os resultados obtidos. A caracterTstica chave, indi_ 
cada, eaquela que sempre foi utilizada para, em conjunto com duas outras 
constantes da Tabela, fornecer o resultado assinalado. 

Varias combinagoes de caracterTsticas de textura fornec^ 
ram urn resultado de 87 acertos, em 91 possTveis, dando uma percentagem de 
classificagao correta de 95,6?^, que foi o maxi mo obtido com caracterTsti_ 
cas deste tipo. 0 nielhor resultado obtido foi de 90 acertos (98,9%), com 
0 conjunto formado pel as caracterTsticas XIDM-M, SUMENT-M e 7-M. Apresen_ 
tam-se, a seguir, as Tabelas de IV. 1 a IV. 11, mencionadas anteriormente. 


TABELA IV.1 

NOMERO DE JANELAS CLASSIFICADAS CORRETAMENTE 
UTILIZANDO-SE CARACTERTSTICAS INDIVIDUAIS 


RACTERlSTICA 

MIN-MAX 

VMP 

CARACTERTSTICA 

MIN-M/ 

ASM-M 

62 

66 

SUMENT-D 

48 

ASM-D 

69 

72 

SUMENT-F 

49 

ASM-F 

68 

65 

ENT-M 

63 

CONT-M 

48 

52 

ENT-D 

41 

CONT-D 

52 

57 

ENT-F 

44 

CONT-F 

48 

54 

DIFVAR-M 

33 

COR-M 

61 

54 

difvar-d 

39 

COR-D 

46 

54 

DIFVAR-F 

41 

COR-F 

42 

47 

DIFENT-M 

59 

VAR-M 

51 1 

42 

DIFENT-D 

53 

VAR-D 

31 

34 

DIFENT-F 

54 

VAR-F 

32 

39 

CORME-M 

53 

XIDM-M 

64 

63 

CORME-D 

56 

XIDM-D 

77 

76 

CORME-F 

53 

XIDM-F 

77 

74 

4-M 

51 

SUMAVE-M 

55 

50 

4-V 

38 

SUMAVE-D 

42 

42 

5-M 

65 

SUMAVE-F 

35 

44 

5-V 

39 

SUMVAR-M 

61 

48 

7-M 

79 

SUMVAR-D 

50 

53 

7-V 

60 

SUMVAR-F 

44 

46 

6-M 

81 

SUMENT-M 

61 

65 

6-V 

54 






NUMERO DE JANELAS CLASSIFICADAS CORRETAMENTE, UTILIZANDO-SE PARES DE CARACTERTSTICAS 
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A“9 

W-9 

A-Z 

w-z 

w-g 

W“l7 

W-3Wa03 

W-iN3JIQ 

w-WAJia 

W-1N3 

W-iN3WnS 

w-yvAwns 
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Q-WQIX 

w-waix 

w-yvA 

w-yoo 

W-iNOO 
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CO 
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VO 
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LO 

VO 

VO 

VO 

LO 

VO 

p^ 

00 

00 

CO 

VO 

00 
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cn 
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LO 


VO 

CO 


00 

p^ 

LO 

VO 

CM 


0 

CM 

LO 
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CO 

00 
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CO 

00 
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00 

CO 

00 

00 

00 

CO 



00 


cn 

VO 

CO 
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cn 
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00 

00 

00 

10 

CO 
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VO 
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LO 
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00 
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VO 
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VO 

0 
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00 

LO 

CO 
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CO 

LO 

CO 

CM 

cn 
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LO 
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LO 
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p^ 

CO 
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0 

, 

r— 

00 

P^ 

0 

CM 

LO 

LO 
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P-H 

VD 
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VO 

VD 
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CO 
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, 
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VO 
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00 
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p^ 
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UO 
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TABELA IV. 3 

CARACTERTSTICA CHAVE: ASM-M 



21 

o 

1 {_ 

2: 

1 

2: 

1 


2Z 

1 






1 

1 

1 

< 

2: 

2: 

LU 

■5T* 






Q 

Q 

A* 

O 

s: 

(jlI 

2: 

1 

q: 

2: 

1 

>• 

1 

s: 

1 



X 

X 

X 


CO 

UJ 

CD 

CD 

} 

r-~ 

1 

t 

CO 


VAR-M 

81 

80 

79 

81 

75 

75 

86 

87 

79 

86 

- Iv 

XIDM-M 


83 

82 

83 

77 

77 

85 

88 

79 

85 


XIDM-D 



77 

82 

78 

78 

85 

86 

82 

85 

>; u 

XIDM-F 




81 

77 

77 

84 

85 

81 

84 

ftp-, 

SUMVAR-M 





76 

76 

81 

88 

77 

86 

V. 

SUMENT-M 






62 

80 

88 

71 

86 

L> 

ENT-M 







•80 

88 

71 

86 


CORME-M 








87 

78 

85 

- ■ 

7-M 









86 

87 



7-V 


85 


. 1 . 


TABELA IV. 4 


CARACTERTSTICA CHAVE: VAR-M 


T'l 

o 

1 1 

s: 

1 

s: 

1 

1 


s 

1 





1 

U— 

uu 

c 

1 

2: 
1 1 1 

s: 

1 

LU 





o 

o 


LU 

s: 

1 

OL 

s: 

1 

>- 

1 

s: 

I 


X 



CO 

ib* 

LU 

o 

1 

1 

1 

CO 

XIDM“M 

84 

83 

79 

81 

81 

75 

84 

78 

86 

\ \ XIDM^'D 


80 

82 

81 

81 

85 

86 

82 

85 

XIDM-F 



80 

80 

80 

84 

85 

80 

84 

i •. 

^ ■ SUMVAR-M 




80 

81 

79 

86 

74 

85 

SUMENT-M 





75 

87 

88 

79 

86 

ENT-M 






87 

88 

79 

86 

1 CORME-M 







83 

81 

85 


7-M 86 86 

h : 

t • 7-V 85 

pi 



TABELA IV. 5 


XIDM-D 

XIDM-F 

SUMVAR 

SUMENT 

ENT-M 

CORME-I 

7-M 

7-V 


CARACTERTSTICA CHAVE: XIDM-M 



82 85 84 84 85 87 84 86 

84 83 83 84 86 83 85 

M 83 83 84 88 79 86 

M 77 87 90 79 87 

86 89 79 86 

1 84 83 85 

87 87 


85 
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TABELA IV. 6 

CARACTERTSTICA CHAVE: XIDM-D 



s: 

1 

S 

1 

U— 


s: 

1 





<c 

e: 

t • f 

s: 

1 

UJ 





s: 

ULJ 

S 
•— > 

1 

H- 

a: 

2: 

1 

>• 

1 

s: 

1 


(/) 

oo 

LU 

o 

1 

1 

1 

XIDM-F 

81 

78 

79 

85 

86 

82 

86 

SUMVAR-M 


82 

83 

85 

87 

83 

86 

SUMENT-M 



80 

87 

87 

83 

87 

ENT-M 




87 

87 

83 

87 

CORME-M 





86 

84 

86 

7-M 






85 

86 

7-V 







85 






' i 






TABELA IV. 7 


CARACTERTSTICA CHAVE: XIDM-F 



s: 

1 


s: 





r~ 
2: 
i 1 1 

s: 

1 

LU 





ULI 

s: 


or 

1 

>• 

I 

s: 




CJ 

1 

1 

r>. 

1 

SUMVAR-M 

81 

82 

85 

86 

81 

85 

SUMENT-M 


79 

86 

86 

82 

86 

ENT-M 



86 

86 

82 

86 

CORME-M 




85 

83 

85 

7-M 





84 

85 

7-V 






84 


TABELA IV. 8 

CARACTERTSTICA CHAVE: SUMVAR-M 



ENT-M 

CORME-M 

7-M 

> 

1 

6-M 

SUMENT-M 

76 

82 

89 

76 

86 

ENT-M 


82 

89 

77 

86 

CORME-M 



85 

77 

84 

7-M 




87 

86 

7-V 





86 


■i! 


I 




i 

I 

1 





J 
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TABELA IV. 9 

CARACTERTSTICA CHAVE: SUMENT-M 


s: 



UJ 

s 

1 

>- 

1 

2: 

1 


o 

1 

1 


ENT-M 

82 

89 

71 

87 

CORME-M 


89 

79 

88 

7-M 



87 

89 

7-V 




85 


TABELA IV. 10 

CARACTERTSTLCA CHAVE: ENT-M 


I I I 

CORME-M 88 79 87 

7-M 87 88 

7-V 85 
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TABELA IV. n 

CARACTERTSTICA CHAVE: CORME-M 


CONJUNTO USADO 

N9 DE ACERTOS 

CORME-M, 7-M, 7-V 

86 


CORME-M, 7-M, 6-M 

82 


CORME-M, 7-V, 6-M 

85 


7-M, 7-V, 6-M 

87 



De unia maneira geral houve um melhor desempenho com a iji 
trodugao de uma terceira caracterTstica, o que pode ser visto confron_ 
tando-se a Tabela IV. 2 com as seguintes. Por exempio, opar ASM^MeVAR-M 
acertou 75 vezes e, com a introdugao de mais uma caracterTstica, obtiv£ 
ram-se resultados maiores ou iguais a 75, o que pode ser visto na la. 
linha da Tabela IV. 3. Isto acontece praticamente com todos os pares de 
caracterTsticas de textura, onde a nova caracterTstica introduzida traz 
bastante informagao sobre a separagao das classes 2 e 3. Nos pares mi£ 
tos de caracterTsticas isto, entretanto, nao e verdade: por exempio: o 
par ASM-M e 7-M que apresentou 88 acertos, com a introdugao de mais uma 
caracterTstica, diminuiu ou manteve este desempenho (Tabela IV. 3, 3a. c£ 
luna da direita para a esquerda). Isto se deve a dois efeitos antagon_i_ 
cos envoi vidos: um, que e a informagao adicional contida na caracterT^ 
tica, e outro, que e a sensibilidade demonstrada pelo MTn-Max, a um au_ 
mento de dimensional idade. No Apendice A. 1.1 sera comentado este seguin 
do efeito. 


De uma maneira geral, pode-se considerar, pelo menos, r^ 
zoaveis os resultados obtidos com trincas dessas caracterTsticas de te]< 
tura. 0 mTnimo de acertos foi 62 {68%), o maximo foi 87 (95,6%), e a mi 
dia foi em torno de 81 acertos (89%). 
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Etnbora aumentando o tempo consumido, o aumento de dimens io 
nalidade do espago de caracteristicas nao afeta o desempenho do classifi^ 
cador,\baseadd na regra de decisao do viEinho mais proximo (VKP). Na 
bela IVJ2 estao alguns resuUados obtidos, quando se utiTizou a dimensao 
do espago igual a 4* Neste caSO, fez-se uso de caracterTsticas que deram 
os mel bores resuTtados, quando usadas em trincas. 

Como se pode observer na label a IV. 12, o con junto de cara£ 
terTsticas de texture ASM-M, VAR-M, XIDM-F e CORME-M obteve 100?$ de cla£ 
sificagao correta, quando o VHP foi usado. Na verdade, este resultado fiao 
se torna muito raro, quando sao utilizadas dimensoes maiores que 6. 

A esta altura, talvez seja interessante comparar alguns r^ 
sultados encontrados, com outros existentes. Weszka et a1 (1976) fazemuma 
analise comparativa de medidas de textura, com fins de classificagao de 
terrenos. Embora as medidas realizadas sobre as matrizes de coocorrencia 
sejam utilizadas, diretamentej como caracterTsticas de textura (havendo, 
portanto, sensibil idade a orientagao), pode-se notar dois fatos: 

1) A entropia (ENT) deu resultados similares aos obtidos com o m£ 
mento angular de 2a. ordem (ASM), quando utilizados individual_ 
mente. No caso deste trabalho, isto tambem ocorreu com relagao 
as medias destes tipos de medidas (label a IV.l). 

2) 0 XI DM deu resultados piores que o ASM e o ENT, quando usados i£ 
dividualmente, ou em conjunto formados por urn mesmo tipo de medj_ 
da, computada para varies distancias e diregoes. Como so se usou 
distancia unitaria, nao se pode confirmer a segunda parte deste 
resultado. Quanto a primeira parte, pode-se ver, pelos resul tados 
apresentados em Weszka et al (1976), que o desempenho do XDM di_ 
minui nitidamente com a distancia utilizada, embora, mesmo com 
distancia unitaria, ele ja seja inferior aos desempenhos do ASM 
e do ENT. No presente trabalho, se se considerar o XIDM-M como 
sendo a caracterTstica mais proxima das medidas diretas sobre as 
matrizes, pode-se ver que os resultados nao diferem muito (os de 
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setnperihos do XIDM-M, ASM-M e ENT-M foram 64, 62 e 63, respecti vamen 
te), einbora o XIDM-M tenha sido uni pouco melhor. Entretanto, o 
XIDM-D e 0 XIDM-F, individualmente, apresentaram os melhores desem 
penhos de todas as caractenstfcas de textura implementadas . Essa 
diferenga talvez seja devido ao fa to de que uma caracterTstica de 
textui^a, do tipo desvio padrao ou faixa de variagao, estamuito mais 
relacionada a direcionalidade, do que, propriamente, a rugosidade 
da textura, que e o que esta, basicamente, sendo medido atraves 
das caracterTsticas usadas em Weszka» 1976. Ou seja, pelo fato do 
XIDM-M ter tido uni desempenho fraco, nao necessariamente o XIDM-D 
ou XIDM-F devem fornecer maus resultados. 

TABELA IV. 12 

CONOUNTO DE CARACTERXSTICAS 


CONJUNTO DE CARACTERTSTICAS 1 

MIN-MAX 

VMP 

ASM-M, VAR-M, XIDM-D, CORME-M 

86 

90 

ASM-M, VAR-M, XIDM-F, CORME-M 

83 

91 

ASM-M, VAR-M, CORME-M, 7-M 

86 

89 

ASM-M, XIDM-M, XIDM-D, SUMVAR-M 

83 

87 

VAR-M, XIDM-M, XIDM-D, SUMVAR-M 

84 

85 

VAR-M, ENT-M, CORME-M, 7-M 

87 

89 

1 XIDM-M, XIDM-D, SUMVAR-M, SUMENT-M 

85 

87 

1 XIDM-M, SUMENT-M, CORME-M, 7-M 

89 

88 

I XIDM-D, SUMENT-M, ENT-M, CORME-M 

87 i 

86 

1 

XIDM-D, ENT-M, CORME-M, 7-M 

I 

86 

90 j 

( 

SUMENT-M, CORME-M, 7-M, 6-M 

88 

85 

SUMENT-M, 7-M, 7-V, 6-M 

87 

87 



4.4 - CONCLUSOES 


De maneira geral podem ser considerados bons os resultados 
obtidos, utilizando”Se apenas caracterTsticas de textura na discHminagao 
de tipos diferentes de terreno. Mesmo tendo sido utilizada uma imagem em 
que nao era tao obv’ia a separagao visual, atraves de textura, de uma das 
classes (tipo C), os Tndices de classificagao correta chegaram a atingir 
varias vezes 95,6%, quando o MTn-Max foi utilizado, e outras vezes cheg^ 
ram a 100% com o uso do VMP, Os tipos de medidas de textura ASM, XIDM, 
SUMENT, ENT e CORME foram os que, em geral, forneceram os melhores resul_ 
tados. 0 CORME, que ainda nao havia sido testado nos trabalhos anteriores, 
forneceu uma boa opgao para distingao entre as classes B e C,asquais, em 
geral, nao foram bem separadas com caracterTsticas de textura individuais 
(Figuras IV. 2 - IV. 6) 

Globalmente, os melhores resultados ficaram por conta dos 
conjuntos mistos de caracterTsticas , ou seja, os que utilizavam tanto c£ 
racterTsticas espectrais como de textura. Esta conclusao esta de pleno 
acordo com o que foi obtido em Haralick e Bosley (1973a). Entretanto, s^ 
ria bom frisar que o uso de conjuntos mistos de caracterTsticas esta re^ 
trito a discrimi nagao entre categorias dentro de uma mesma imagem, uma vez 
que as caracterTsticas espectrais de uma mesma classe podem variar sensi^ 
velmente, de acordo com as condigoes sob as quais a imagem foi obtida.De^ 
te fato e que decorre a grande vantagem das caracterTsticas texturais so 
bre as espectrais, uma vez que antes do processamento e realizada a equ^ 
lizagao do histograma, que as torna invari antes a condigoes diferentes, que 
porventura existam em imagens distintas. 
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PROCEDIMENTQS DE CLASS IFICACAO 

No programa CATEG, desenvol vido para realizar a etapa de 
classificagao propriamente dita, podem ser utilizados dois classificado 
res: o MTn-Max e o do vizinho mais proximo, 

A.l - 0 CLASSIFICADOR MTN-MAX 


Seja uma amostra designada por urn vetor de N caracterTsti^ 
cas X = (xi, X 2 •••» x^). Este classificador utiliza uma regra de decisao 
que e de maxima verossimilhanga sob a hipotese de que as caracterTsticas 
sao independentes e possuem distribuigao uniforme. Se isto for admitido, 
tem-se que a fungao densidade de probabil idade para a categoria k e: 

f{Xlk) ^ , para todo X tal que 


‘’nki Xn-Sk- 


onde, nesta expressao, a^j^ e b^j^ definem os valores maximo e mTnimodadi^ 
tribulgao uniforme da categoria k, para a componente n. 


Assim, uma amostra (xi,X 2 ,.,., x^) e atribuTda a categoria 


k, se e so se: 


1) ^nl ®nk’ N 


2) ^ , para todo j, tal que 


^^nk " *^nk^ 


(^nj ■ 


^nj - ^n - ®nj’ ^ =1,2,..., N 
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Se nao existir k, tal que £ a^l^, n=l,2, N, 

entao, X e atribuTdo a categoria k, se esta for a categoria que minimize 
0 termo 




onde K e 0 numero de categorias e o (*) foi colocado para indicar que nao 

se incluem no somatorio, os termos cujos n'S sao tais que b • < x„ < a .. 

nj — n — nj 

A estimativa do intervalo b^j^ - a^j^ da distribuigao unifor. 
me da categoria k, para a componente n, e feita da seguinte forma: Sejam 

Ank 6 os valores maximo e mTnimo da componente n, para todas as amo£ 
tras designadas na categoria k, e Mj^, o numero dessas amostras. Entao: 


^nk ' ®nk 


\k - ^nk 
M, -1 


%k = ''nk ^ 


*^nk - '^nk 

Mk-1 


Ve-se que o intervalo da distribuigao uni forme e urn pouco 
maior que o encontrado na fase de treinamento, isto e: 

M. _ , 

^nk ■ ^k " “ ^k^ • 

Mk - I 

Nota-se que e bastante simples o funcionamento do Min-Max. 
Na fase de treinamento, e determinado o N-paralelepTpedo correspondente a 
cada uma das categorias existentes. Quando, na fase de classificagao, uma 
amostra e apresentada ao classificador, este verifica em qual paralelep^ 
pedo el a esta contida. Podem entao ocorrer 3 casos: 
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1) 0 paralelepTpedo e unico; entao aamostrae classificada como se]i 
do da categoria correspondente. 

2) 0 paralelepTpedo nao e unico: entao a classe correspondente ao 
nor deles e atribuTda, pois, de acordo com a hipotese feita, el a 
e a mais provavel . 

3) A amostra nao pertence a nenhum paralelepTpedo. Neste caso, pr£ 
cura-se aquele paralelepTpedo que esta mais perto da amostra, e 
a classe desse e atribuTda a amostra. 

0 terceiro caso envoi ve calculo de distancias euclideanas, 
no espago de caracterTsticas . Como, certamente, as caracterTsticas nao tern 
a mesma ordem de grandeza, e necessaria uma normalizagao, demodoase atr£ 
buir pesos iguais para todas el as. Foi feito urn mapeamento linear em cada 
dimensao do espago de caracterTsticas, i ndependentemente da classe das 
amostras, de modo que os valores assumidos variassem entre 0 e 1. 

Quando da utilizagao do MTn-Max, o conjunto total de amo£ 
tras foi dividido em duas partes iguais: 91 amostras para treinamento, ou 
determinagao dos parametros das distribuigoes pelo cl assi fi cador , e 91 
amostras para teste. Foi verificado que, quando era aumentada a dimensao 
do espago de caracterTsticas, crescia o numero de ocorrencias de classifi_ 
cagoes, pelo terceiro criterio citado, ou seja, mais amostras nao perteji 
ciam a nenhum dos paralelepTpedos existentes. Isto talvez seja consequiji 
cia de uma estimativa do interval o da distribuigao uniforme de cada cla£ 
se, para cada caracterTstica, baseada em urn numero insuficiente de amo£ 
tras. Como no caso deste trabalho as classes nao eram separadas, o uso d£ 
quele criterio era o que mais resultava em erros de classificagao. 

Para simplificar uma analise, considerar-se-a o caso em 
que 0 intervalo da distribuigao, uniforme em cada componente, e estimado 
como sendo a diferenga entre os valores maximo e mTnimo das amostras para 
esta componente. Supondo-se que se faga tal estimativa, baseados em duas 
amostras apenas, e X 2 , e que as caracterTsticas sejam realmente ind£ 
pendentes e uniformemente distribuTdas entre 0 e 1, pode-se, entao, calc£ 


lar a probabilidade de que, pelo menos, uma fraqao w das amostras venha 
cair fora do intervalo estimado, no caso de dimensao igual a 1 (FiguraA.l), 



Fig. A.l - Fungoes densidade de probabilidade. 


A largura do intervalo estimado e t = |xi - X 2 I, e a pro 
babilidade que se procura e (Figura A.l). 

P Cl - |Xi - X 2 I > W] = P Cl Xi - X 2 I < 1 -w] 

Calculando-se este valor, atraves da fungao densidade de 
probabilidade de variavel aleatoria |Xi - X 2 I, apresentada na Figura A.l, 
tem-se: 

fl -w 

P CiXi - X 2 I 1 1 -w] = (2-2y) dy=l -w^ 

0 


Calcule-se agoy’a a probabil idade de que, pqlo menos, uma 
fragao y das amostras caiam fora do intervale estimado, no case de dimen 
sao igual a 2 (Figura A. 2). 


Fig. A. 2 - Intervaios bidimensionais . 


A partir da Figura A. 2, e levar lo-se em conta o resultado 
do anterior, tem-se que: 

P C 1 -IXia - X 2 b I la. 1 -|Xib - X 2 bl 1 bj = (1 -a2) (1 -b2) 

e a probabil'^dade de que, pelo menos, uma fragao 1 -(1-a) (1-b) das amos^ 
tras caiam fora da regiao estimada. 


Fixando-se uma certa probabil idade, e comparando-se as fra 
goes mTnimas esperadas de amostras, caindo fora do interval o estimado p^ 
ra os casos de dimensao 1 e 2, pode-se ver que e realmente esperadoumnu 
mere maior de amostras fora da regiao estimada, para o caso de dimensao 
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' 2. Por exemplo, com uina probabilidade 90$^, espera-se que, pelo menos, 31 ,65» 
das amosti^as caiam fora para o caso de dimensao 1, ao passo que, com dj_ 
mensao 2*, esta percentagem mTnima se eleva a 40,25^, 

Para uma estimativa baseada eni urn numero muito grande de 
amostras, a situagao nao e tao simples assim, Porem, e razoavel que ainda 
neste caso, auniente o numero de amostras caindo fora do intervalo, para 
urn aumento de dimensional idade, uma vez que, para uma amostra cair dentro 
de uma regiao a N dimensoes, ela precisa estar dentro de cada urn dos N iji 
terval os corres ponden tes . 

Como se pode verificar na descrigao do MTn-Max, a largura 
do intervalo estimado da distribuigao uni forme, para cada componente, nao 
e exatamente a encontrada atraves das amostras de treino, pois existe urn 
fator mul tipi icativo maior que a unidade, e que e fungao do numero dessas 
amostras. Entretanto, ainda assim foi verificado urn cresc^mento do numero 
de amostras em teste caindo fora do paralelepTpedo estimado, com o aumen^ 
to de dimensional idade. 

A. 2 - REGRA DE DECISAO DO VIZINHO MAIS PROXIMO 


Seja urn conjunto de amostras, com classificagao conhecida, 
{Sx, S 2 S|^}, e X a amostra sendo classificada. A regra de classifica 
gao do vizinho mais proximo consiste era atribuir a X a classe do seu vi_ 
zinho mais proximo, s*, onde s* e tal que: 

D (s*, X) = mTn {D(s^,X), n = 1 ,2, . . . ,N} 

Nesta expressao, D pode ser qualquer niedida de distancia. 
No caso deste trabalho, a distancia utilizada foi a propria euclideana. 





De tii0do a se atribuir a mesma importancia a todas ascarac 
tensticas, deve-se realizar uma normal izagao, de maneira que todas ass£ 
mam valores dentro de urn mesmo interval o. 

Este classificador, embora mais elementar que o Min-Max, 
apresenta a desvantagem de ter uma velocidade de classificagao menor e 
mais dependente da dimensao do espago de caracterTsticas. 
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C PROGRAMA PAHA aRTb'NCAO SJSTKF.AnCA DE GAHACThR ISTICAS 
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0001 
0002 
0o03 
• 0004 
0005 
000b 

0007 
000 8 

0009 

0010 ■ 

0011 5 

0012 

0013 .10 

0014 rlOSS 

0015 vOOO 

0016 
0017 
001 8 
001 ') 
if </ K i. 

0023 

0024 

0025 9001 
.0026 

0027 

0028 9002 

0029 

0030 
00 il 

0032 

0033 

0034 9010 

0035 

0036 . 
0038 

0040 12 

0041 . 1 1 

0042 

0043 

0044 20 

0045 
004b 
0047 
004(3 
0049 

0051 30 

0052' . 

005.3 40 

0054 


LOGICAL?! lvC74) , IZ , TA(3I,K( 256 ) , &UFF( 5 1 2 , 2 ) , RUF ( 5 1 2 ) 

BEAL»4 VPW(12)#VPiv(l2),XVAP(4),XMEDlA(4) 

PEAL ♦A VQW(12) ,FEDEPU2) ,FEMED(12) ,FERAG(12) , VSW(12) ,VWC12) 
ILTEGEP LXC3,5) ,CX(2,4) ,EC(2,2) ,ICC2,2) ,CC(2,2) 

IhlEGP^R SC(2,2) ,AR(256,4) ,l.r*(2) ,CHK,L>ISl,Lt(vA(4) ,NURHC4) 
C0M:-’'0N P( 32/ 32) ,ALFA(32) ,EsE'U( 64) ,GAF.A(32) 
EOUlVALERCECMi.lZ),(bUF( 1) / RUFK1,1 )) 

EQUIVALENCE ( CC C 2 , 2) , A’AP ) , ( CC ( 1 , 2 ) , lAR) , ( CC ( 1 , 1) ,KAR) 

CALL OUTPUTC27/12) 

1-.P1.TEC6/5) 

FnRt>)ATC13X, ' ♦ + ? EXTHACAO L'E CARACTERISTICAS i , // ) 

WRITE (6, 10) 

fORMAj'dX., ' A IMAGEM DEVE TER SIDO CAPBEGADA USANDO A MLSMA 
C ESC A LA DO ERTS') 

WRITE l 6/9fi00) 

FuRJ^AvC < sprtiMfc Ih-A VEZ QUE USA 0 PROGRAMA { E/'N ) ?> ' ) 

CALL UU1PUTC7) 

REA0C6,20)W 
CALL FPOf<(T(W,74) 

1 r w I 1 . EO . " 1 ! 6 } GO TO 90C ’ 

It I wi j ) • IMS, , i j I jU lU b9l*y 
CALL FD6SEK3) 

GO TO 9010 
lwED=l 

CALL FDBSETC3/ 'OLD' ) . 

WRiTE(6/9002) 

FO.RKATC ' SQUANTAS JANELAS FOBAM PROCESSAUAS?> ' ) 

CALL OUTPUT! 7) 

READCb,20)W 
CALL FB0;')T(V^/74) 

L = 0 

CALL 1NTFFCIj,W,74,.'4PEC) 

CALL ASSIGNC3, 'OTOzTAE.OAl ' ) 

DEFINE FILE 3 C 200 , 1 00 , U , NREC ) 

IF ( IwED.EO. 1 )LPEC=MPEC+1 ’ 
iKClWEU.NE.n 1'1REC = 1 
WRITE (6,11) 

FORMAT! 'SDISTANCIAC 1-33 t CANAL! 1-4) A USaR?>') 

CALL. OU'fPUT(7) 

READ!6,20)W 
FG.''M,.A'H74A1 ) 

C ALL FRuN’i ( w , 7 43 
L=0 

CALL INIFF CL,W,74,01ST) 

CAf.L JN!)FF!L,lv,74,CHK) 

lF*(DIST.Gr.3.0R.CHM.GT.‘4.0R.01Sl.LT.1.0P.CHN.LT. l)GO. TO 12 
DO 50 J = 1 ,2 
WRITE! 6/ 40) J, J 

forma 1 ! ' S b at a ex ' » I l . ' » El ' • X 1 f ' ; ' f 2X ) 

CALL OUTPUT! 7 ) 


FORTRAfJ 

CORBsOaK 


IV V01B.-02 

, UJC=t200,200J. 
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CQbbTA , L.P ; / LI : I =0'1 1 : COLBTA . 


0055 

* 

REAO(b,20)W 

0056 


CALL FRONT(W*74) 

0057 

‘A 

IF(Vv(l).EQ.''130)GO lU 2000 

0059 

i 

L=0 

0060 


DO SO K=1 ,2 

0061 

50 

CALL lf!TFF(L,W.74,EClK,J)) 

0062 


IF(EC(l,2) .LE.ECn , U .OR.LC(2,2) ,LE.E;C(2, n )GO 

0064 

60 

DO 80 J=l,2 

0065 


WRn'eib,70)d,J 

0066 

* 70 

FORMATC 'S6ATA IX ' , 1 1 , MY ' r 1 1 , ' : ' » 2X) 

006? 


CALL OUT pm (7) 

0068 


READ(6,20)W 

0069 


CALL FPOWr'^W,74) 

0070 


rF(W(l).EO.''130)GO TO 2000 

0072 


L=0 

0073 


DO 80 K=l,2 

0074 

80 

CALL TMTFFCL,U«,74, 1C(K,J) ) 

0075 


IFnCn,2J .LE.lCCl , 1 ) .OP.XCC2,2) .LE.ICC2, l))GO 


'Z tSCOr.KA E DETERrilNACAO DOS PARAMETHOS DA JAKbLA 

0077 

* B9 

CALL JAWe.LA(EC,lC,CC,SC,JOb,IDX, IDY) 

)G78 


iFCJ0a.E3.0)G0 TO 2000 

GO 80 

27 

WRnE(6,‘50) 

0-00 1 

•JO 

FORMAl ( ' SQUAL A CATEGORIA DA JANtLA bSCOLrill.>A?> 



CALL 0UTPUT(7) 

008 3 


R£AD(6,20)W 



fAf.i. Htnurri-. ,7 41 

U V J 


L-'j 

00^6 


CALL IK‘TFF(L,w,74, ICLASS) 

00 87 


IFCICLASS.LT. 1 )GO TO 27 

0089 


CALL ASSIGNCS, 'LP; ' ) 


C OBTENCAQ DO H^STOGHA^'A•• E CAPACTLPISTICAS ESPECTRAIS 

0090 


DO 120 J=1 ,4 

0091 

* 

LXU ,J)=0 

0092 


LX(2,J)=0 

0093 

120 

•LX(3rJ) = l 

0094 


. LX(l,b)=-l 

0095 

* 

CALL IWHL(LX) 

0096 


DO 150 J=1 ,255,2 

0097 


DO 140 1=1,4 

009B 


DO 140 K=l,3 

0099 

140 

LX(K,l) = LX(K,.l)+2 

OiOO 


CALL IrtHC(CX) 

0101 


CALL IivHLCLX) 

0102 


DO 145 M = l,4' 

0103 


DO 145 K=0,1 

0104 

145 

AP(J + K,MJ=CX(S=H ,M) 

0105 

150 

CO.’mNUE 

0106 


CALL IPTOTCLENA) 

0107 


GALL IRTOT(LEhA) 

0108 


XAREA=R101 (LENA,2) 

0109 


DO 170 M=l,4 

0110 


XMEDI A(M)=0.0 

0111 


XVAR(M1=0.0 

0112 

' 

DO loO J=2,256 

A 1 •> 
U 1 V ^ 


ACUM=RUS( ARC J,M) ) 

0114 


XMEDIA ( !•:) -XMEDI.A ( M ) + ( J- 1 ) ACUM 


ii 

II ^ 

<4 a 
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CQRB5 

oaK» 

UlC«[2|iO;y200J CQUE1'A,LP:/L1 Sl-DTl J.COEETA. 

Oils 

160 

XVAK(M)=XVAIUM)+IJ’*1 ) 

0H6 


XKEOl A ( M ) aXMEDl A ( H) /XAREA 

0117 


XVAH(«) = XVAK(K)/XAREA-Xf-'EU1A(M)*XKED1ACN) 

01 la 

170 

COMTINUE 

* 

c fc.Olj'Al.lZACAO DU CAK’AL BSCULHIDU EM 32 NlVblS PARh D CA^AiJ b 

ouv 


ARbA=0.0 

0120 


DO 173 Jsl,256 

0121 

173 

AHEAsAkEA + RUS(AR(0,CH.N) ) 

0 122 


y-RlTE(6,17b)CriN 

0123 

* 175 

FORMA'K lOX, ' iraCIO DA EOUAEIZACAO DO CANAL', 12) 

0124 


ACOM=RUS(AR(l,ChN)) 

0125 


MlaACon+ 32/AREA 

0126 


U (NJ .EO.O)GO TO 176 

012 8 


TnrA=AREA-Fp!l/32 

0129 


If (ACUM.Gf .TOTA+-AREA/64)GO TO 176 

0131 


Nl=in-1 

01 32 

176 

KsNX 

013 3 


Nl=Nl+l 

0134 

* 

TABLE! UalZ 

013^ 


DU lyO 0=2,256 

n 1 Xh 


Ktni=FHlS(AR(.l.CHN) ) 

0137 


. ACij'i;=ACu*'5T&ADi 

0139 


(.1 1 rzpQUUt 3 2 / ARE A 

0139' 


IFCM.EQ.KIGO TO 180 

0 1 4 1 


rriT’A=ARFfl*Kl /32 

■»M T K » - rt > - » > » < 

»./ 1 4 2 
0143 

* 

1 A n — c, w n u X 

IF( (rEIA + ACUHJ'/2.GX.'JUTA)G0 TO 1/8 

0145 

• 

Nl=NI-l 

0146 

178 

K=)a • 

CM 4 7 

180 

MI=NI+1 . • 

0148 


IF(N1,GT.32) NI=32' 

0150 

190 

XABLE(J)=IZ 

0151 


LIVV=CC(2, l)-3 

0152 


LSVV=CC(2,2)iF3 

0153 


LIVH=CC(1 , l)-3 

0154 


LSVH=CC(1 ,2)+3 

0155 


DO -220 L=LIVV,LSVV 

0156 


CALL IRV(CHN,L,BgF) 

0157 


CALL aAIT 

0158 


DO 210 J=LIVH,LSVH 

0159 


k=IbyfE( J-1 ,BUF)+1 

01 bO 

210 

BDflJ)=TABLE(K) 

Olol 


CALL I/<V(5,L,faUF) 

01b2 


CALL WAIT 

0163 

220 

CONTINUE 


C OBTEfJCAO DAS i.lATRIZES OE COOCORRENCIA 

0164 


N.N=CCCl, n+DlST 

0 1 65 


(.•UR-4 ( n=2» ( lOX-DlST) » IDX 

0166 


NCi.Rr.U2) = 2+( 1DY-DIST)*IDX 

0167 


M 0 P M C 3 ) = 2 » U 0 X - D 1 S X ) * U n Y - D ). S T ) 

0 3 o8 


NORf'U43=fyORH(3) 

01 b9 


MODOaj 

0170 


NBn)*2 

0171 


Wb(2)=l 

0172 

240 

DO 255 J=l,32 

0173 


DO 255 1=1,32 


FORTRAN IV VOl‘i-02 
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0174 

0175 
01 7b 
0177 
•0 17 y 
0179 
01 BO 
0181 
0182 
018 3 
01 8 4 
0185 
0106 

0187 

0188 
0190 
0192 

19 3 

0194 

0195 

0196 

0197 

0198 
0200 
0202 


02U4 

0205 

0206 
0207 
020 8 
0210 
0212 

0213 

0214 

0215 
021b 

0217 

0218 
0220 
0 222 

0223 

0224 
0 2 25 
0226 
022 7 


0228 

0229 

0230 

0231 

0232 
02 3 3 
02 34 
0 2 3 5 


255 PCl.,J)-0.0 

DU 401 JKsl.DiST 
Ib=l 

Jli = CC(2; 1 )+JK-l 
JM=JtifDlST 

CALL lHV{5,JL,BUFF(l,Ib5) 

CALL WAIT 

00 400 KsJMrMAI^DlSr 

GALL lHVC5,K,bUF F( l,Nb(lb))) 

CALL 

GO 10(300,305,310^ 315 )FODC 

300 DO 301 .N = N'M,IA1^ 

I=n3Yl£(N-DIST, 6UFP( 1 , IB) ) 

JrltiYlE(M,&UFF Cl , 18) ) 

JF(1.GT.32) I=TABLEC1) 

IF(J.G1\32) J=1ABLECJ) 

PC 

301 FCJ,I)=PCJ,1)+1. 

GO TO 40Ct. 

305 DO 302 N;^KAH,IAR 

l-lbYXh.(!\',BUFFCl , 18) ) 

J = 1 8 \ T £ C , 5 0 F F C 1 , N B Cl B ) ) ) 

IFM.Gr.l?) r = TABLB(l) 

IFCJ.GT.32) J=TArtb£CJ) 

PCI*. J)=P( 1, J3XU 

1 ‘\ -4 n t • ! . — r. 4* 1 T •. * 1 

-4 4* 44. 4V444r*^ *\4^f«4«*4<4 » 

OU lu MOO I 

310 DO 303 K=NN,IAB 

J = lb'tTL-C(v-D,iSl ,.B’JFF(1 ,F'Ba&))) 

1 = 1BYTE.(N,BUFF(1,I3)) 
lFCl.yT.32) I=.TAbLE(I) 

1FCJ.GI.32) J=TAbL£(J) 

P(1,J)=P(1,J)+1. 

303 PCJ,1)=P(J,I)+1. 

GO TO 400 

315 DO 304 N=NN,1AB 

2=1BVTFCW~D1ST,BUFF(1 ,1B)) 

J=ibYTE(N,BUFF( 1 , NBCIB) ) ) 

IKI. G1.32) I=TABLE(1) 

IFCJ.GT.32) J=TABLECJ) 

PC1,J)=P(1,J)+1. 

304 PCJ,l)=PCJrl)+l. 

■4 00 lB=.‘<aCJB) 

401 COMTINUE 

WRllt Cft,901)wODO 

901 FOkMATClX, 'COMPLETA A OfaTENCAO DA KA’I P 12 ' , 12 ) 

C D0R.V4AL1ZACA0 DAS MATP12PJS 
Du 4 05 J- 1 f 3 2 
DO 405 1=1,32 

405 PCI, J)=PCi.J)/NORH(’4CiDO) 

C UBTtkCAC) DE CA.HACTFHI5T1CAS TE5CTDRA1S A PARTIP DAS MATRJ2FS 
CALL, FtXTCvv.) 

DO 410 1=1,12 

GO TOC406,407,40B,409)FODO 
4 0 0 V p rt C i ) = V A' 1 1 ) 

Go TO 410 


fhl 26-t\iUV“7& 20;19:2j .^AGE’ fU/S 

COGF'iA, LF't/Ll : 1=U'I J tCOLETA. 


FC)Pi'PA4 IV 

com- 

09K, 

0236 

107 

0237 


0238 

406 

0239 


0240 

409 

0 24 1 

4 10 

0242 


02 43 

900 

0 24 4 


0245 


0 24 0 

420 

0247 


0249 


0249 


0250 

. 

0251 


0 25.2 


0253 

450 

0254 

‘ 

025 S 

5 00 

u 2 5 6 


A '1 c ^ 
U 4. ^ r 


025- 


0/59 

2 yj u u 

0260 


0261 

404 

0262 


0263 


0264 



V0li5-02 
UlC=t200,200) 

a 

vd»nn=VA’(i) 

GO TO 410 
VRW(i )=V4(X) 

GO TO 410 
VSrt(l)=VW(I) 

CONTINUE 

>vRlTh.(t.,900)MOOO 
FOtUlAfClX. 'flM DA EXTP.ACAO DE CARAC TERIK’J. IC AS DA KATRIZ',121 
MU1)0=/.UD0+1 
IF ( mODO’4 ) 240 , 240,420 
DU 4S0 1=1 ,12 

FEMED ( I ) = C VP7; C 1 ) + V'JN ( I) +VRW ( U + VSW (1)) /4 
2=FEMEDCI) 

IAB=1VPW(I J-Z)»*2+CV0w(I)-Z)'f + 2+(VRWa)-Z)T + 2+(VStv(l)-Z)t=i=2 
FEDEPIDnSOHTCTAB) 

AMAX = Ai«)AX1 (VPwCI) , VOW ( I ) , VHW (I ) , VSwCU ) 

A^•.lfJ=A^‘^'MltVPwCl) ,VQw(I) ,VPwCl) , MS', HD) 

F’ERAGU1=AMAX-AMIN 

ftplTECbrbOO) C( CF EPSDCl) , F E'DEP ( 1 1 , FERAG ( I) ) , 1 = 1 , 12) , ( (XFEDIAC I) 
CXVARU ) ) , 1=1 ,4) ,SC,1CDASS) 

FQFJ'A n 1 X , 1 2 { 3E1 7 . B , / ) . 4 ( 2EI 7 , 2 , / ) , 5 16 , /// ) 

CALL CtiUSt. (5) 

V’R ITE ( 3 ‘ WREC ) ( ( C FEMED ( I ) , FLDEF ( I ) , FER AGC I ) ) , 1 =1 , i 2 ) , ( ( XMEDIA ( I 
CXVARCI ) ) , 1=1 ,4) ,SC , ICLASS) 

GO rn 69 

!i| w ij s: •• t 

ivt<lTEC6,404)NRD 

POPMATC IX, ‘TOTAL Dr: JANELAS PPOCESSADAS= ‘ , I 4 ) 

CALL CLOSE(i) 

STOP 

END 



WiKJJNAL PAGE TS 
01'^ POOR QUALITY 


IV VOlb-^02 FHJ -26-KOV-76 20:21 :26 PUCE 001 

CORH:«03Kr U1C=[200, 200J * ^ COLKI A , tJ^ : /LI : 1 1 : CuLb • 


0001 
0 0 0 2 
0U03 

0004 82 

0005 

0006 BS 
0 00 7 

0008 83 

000 9 , 

0010 90 

0011 

0012 92 

0013 

0014 20 

0015 

0017 

0018 
0 0 i 9 
0020 
0021 
oov? 

0023 

0 0 2 *t 9 3 

0025 

0 02 6 

t ) 0 2 1 

0028 

0 0 3 0 

.0 0 31' 

00 3 2 
0 0 3 i 
Op34 

003H 95 

0036 100 

003 7 

0 0 3 H 110’ 
00.^9 115 

C040 
0041 . 

00 42 
0 0 4 3 

0 0 4 4 

0 0 4 o 
0048 
0050 


SUBROUl Ifv£ JANELACEC, JC,DD,5C, JOB, lux, IDY) 

1 0 ILCiLR t!.C ( 2 , 2 ) , IC f 2 , 2 ) , DD ( 2 , 2 ) , SC C 2 , 2 ) , CUHSQF C 5 0 
LLGICALM i\C74) 

CALL UU1PUIC27 , 121 
wRnfc:cb,85}H:c 

FORMAT ( ’ EX1 = ‘ ,i6,3X, 'fc:Vl = Mb,10X, 'F.X2=' ,Io,3X, »EY2-M6) 

• HITE (6,68) IC 

FORMAT ( ' IX) ,16, 3X, ' IY1 = M6, lOX, ’ 1X2= * , 1 b , 3 X , ' J V 2= ' 6 ) 

.VKlTECb, 90) 

FUR'^AK/, U, 'LSCOLriA DA JA(vELA:CR P/ LLR COOHDLNADAS EHTS',/) 
riUMi^=i 

CALL OUrpUT(7) 

REAU(6,20)W 
FUkMA r ( 7 4 A 1 } 

IP^(vvU) .t(^,”130)GO TO no 
CALL iRKCCUKSOPJ 
' (1,1 )=CUPSDP(2)-2^CIjRvSQP(3)-l 
UU( 1 ,2 J=CURSCjR(2)+?4CUFSOP(3) 

DJC2 , 1 ) = CilRSaR( 4) -24CUHSOPC5) 

DO ( 2 , 2 ) =C UHSOH C 4 ) + 2 ^ CUBSOk f 5 ) + 1 
DO 9 3 n 2 
Du 93 K=n*2 

SC C K , J ) =DD C X , J ) - IC ( K , 1 ; T L£ ( K, n 
J0X = DD(l,2)-DDCin )+l 
TDy-D*:^\2,2;-DDC2n)kl 

;*• M {XI IJ V M <{- I 

IFCNUMB. LE.13)G0 TO 95 
CALL UUrPUT(27,12) 

.vHilLC6,65)EC 
wHITE(6,88)1C 
;vHITEC6,90) • ’ 

LUM,B = 1 . . 

wklTE(6, 100)SC,IDX,IDy 

FORM At C » Xl=ni4,2X, 'YI = M4,3X, »X2=' ,14,2X, *Y2=n 
CI4,2X,'DX=M3,2X, 'DY=' ,13) 

GO 10 92 
vsHnE(6,115) 

FOKMATC ’ S5FGIJIF ,CICLAR OU PARAP? CS/C/P ) > V, 2X) 

CALL 0U1PUT(7) 

RFAD(6,20)W 
CALL F PO^*TCvV,74) 

JUtt=l 

ir ( a'u ) .Fu.“l03)GU TO 82 

IFCaCI } .rjf:./'123.ANr),W(l) ,LF.<’12G)‘GO TO 110 

lHkv(l).EO*n20) J0a=0 

RETURN 



.FQRTRAIJ IV 

VOlB-02 • FPl 

COKK = 

08 K, 

UlC=t200,2003 

j 

jooui 


SUBROUTINE AUX 

|0002 


COMMON PC32,323 , ALFAC32D 

!000 3 


DO b 1=1,32 

0004 


ALFAU) = 0.0 

•0005 

• 

DO 5 J=l,32 

0006 

5 

ALFA(l) = ALFA(i7 + PCJ,l ) 

U00 7 


RETURN 

0008 


END 


-7t7 J9: li :30 • PAGE 001 

SUBkOTS , LP :/Ll : 1 =DT1 : SUBPOTvS. 


FORTRAN IV VOia-02 

CORE=OBK# UlC=t200,200] 


PHI 26-BOV-76 19 :11 ; 54 PAGE 001 

SUBROTS , LP : / LI: : 1 =0T 1 ; SUBRQTS . 


0001 

SUBROUTINE ASM CPI) 

0002 

COMMON P(32,32) , ALFA C 32) 

0003 

FI =0.0 

0004 

DO 10 1=1., 31 

\r\f\3x? 

Fisf i+pu,i)*Pti^n 

OQi)b 

K=r + i 

0007 

. DO 10 .J=K,32 

0008 10 

F1=F1+2*P(1,J)*PCI,J3 

000 9" 

Fl=FllP(32,32)+PC32,32) 

0 n 1 0 

RETURN 

0 C 1 A 

K u 


FORTRAN IV V01B-02 

CQkE=08K, UlC=t200,2003 


FRI. 26-HOV-76 19:12:18 PAGE 001 

SUBR01S,LP;/L1 : 1=DT1 :SUbPOTS. 


0001 


SUBRDUTINE AUXl 

0002 


COMMON PC 32, 32) ,ALFAC 32 ) 

0003 


DO 35 1=2,33 

0004 


BETA ( 1) =0 . 0 

0005 


DO 25 M=l,32 

000 D 


N = I-M 

0007 


1FCN-M)35,30,20 

C009 

20 

HETACI )=Bt;TACn+2?PCM,N) 

0009 

25 

CONTINUE 

0010 

30 

BKrA(l)=BETA(I)+PCM,\) 

0011 

35 

CONTINUE 

0012 


DU 55 1=34,64 

0013 


8ETA(1)=0.0 

0014 


K=I-32 

0015 


DO 45 ^■«=K,32 

0016 


N=I-M 

00 17 


1F(N-m)55,50,40 

0018 

40 

BElACl)=HElA(l)+2i‘P(M,N) 

0019 

45 

CONTINUE 

0020 

50 

BErACI)=6ETACI)+P(M,N) 

0021 

55 

CONTINUE 

0022’ 


RETURN 

002 3 


END 


original 

OF POOR OUAOUX 


FORTRAN IV V018-02 HU'26-NUV-76 l»;r2!4o PAGE 001 

CORE = 08K, U1C= [2013/1:2001 . RllUROTS , LP : /LI : 1 =DT 1 : SUUROTS , 


OOOl' 

0002 
000 3 

0004 

0005 U) 

0006 
0007 
OGOd 
0 0 00 

0010 20 
OUll 100 
0012 
0013 


SUbRUUTlME AUXO 

COMMON P ( 3 2 , 3 2 ) , ALFA (32), bKTA ( b 4 V, G A M A ( 3 2 ) 
GAi>;A( I )=0.0 
on 10 M=t,32 
GriMAd )=GAMA(l)tH(M,M) 

00 100 K=2,32 
L=32-K+l 
GAMA (K ) =0.0 
DO 20 M=1,L 

G A M A ( K ) = G A M A C K > + 2 ♦ P ( F , K + K - 1 ) 

cooni'vue 

RbTfcJRN 

END 


FORTRAN IV V018-02 HRI 26-N0V-76 10:13:12 PAGE OCl 

CfJpE=OrtK, UiC=[200, 20'0) bilBfil; Xb , l.p: / L 1 : 1 =D T 1 : SUBRO'IS. 


0001 ’ 

0002 

0003 

0004 

0005 10 
0 0 C’ 6 


SOEhuDi I K'E C0r.'i(r2) 

COMMON P(32,32) , ALFAC32) ,BETA(64) , GAMA (32) 
F2=0.0 

DO 10 M=],31 

r X ~ r ? T *•» 'f V, ■>r {.i A m A ( f* -t* 1 ^ 

p r r n D r , 


0007 END 


FORTRAN IV • V018-02 FRI 26-NOV-76 19513:37 PAGE 001 

CCRE = OHK, UIC=[200,200) . • SU8R0TS , LP 5 /LI : 1 =L'l 1 : SUDRUl S . 


0001 


SUBRDUTl'ME C0RCU,F3) 

0002 


common PC32,32J ,ALFA(32) , BFX A C 64 ) , GAMA ( 3 2 3 

0003 


U = 0.0 

0004 


V = 0,0 

0005 


w = 0.0 

0006 


DO 20 1=1 r 32 

000 7 


S0MA=0.u 

OOOR, 


DO 10 J=1 , 32 

0 009 

10 

SOMA=SOMA+J^P(i; J) 

00 10 


U=U+1^ALFA( I) 

0011 


V = v^-f MH^AOFACI) 

0012 

20 

W=vf SOMALI 

0013 


FuK=U*U 

0014 


V=V-FOR 

0015 


F3=Ca-F0R)/V 

0016 


RtrORW 

0017 


EK'O 


IV VOlh-02 

Ci''i?F. = 03K, UIC = [200,2ti0)- 


FFI iti-.'.OV-Tb I9:l-i:03 Ool 

SJbf-’.UJ'vS . L^l/Lltl =!'T1 ; bUPFoi {> . 


0001 

SUUPlUTINE VAR(U,F4) 

■ .•> 

000 2 

GO-t.-lOR PC32, 32) , ALFAC 32) ,«tTA(64) ^(;A.1A(32) 


0003 

F4 = 0.0 


UOO‘1 

00 10 1=1 , 32 


Ooyi, 

H = l-U 

. V 

gooo 10 

F^»=K4 + H + B + ALFAU) 


0007 

. Rfcl'URfv 


0008 

END 



FORTRAN 

IV VOlB-02 FBI 26-NQV-76 19:14:27 PAGE 001 


corh:= 

OBK 

, UIC-C200,200) SUBRQ1S,LH;/LI: 1=DT1 :Sb&P0T5. 


0Qv)l 


SUtfRQin 14c: XIDMCF5) 

... 

0002 


Curi'.ON P(32,32) ,ALFA(32) .BFTAfbA) ,GAMA(A?) 

ii’ 

0 00 


F 5 = 0 , 0 

"1 

0004 


DO 2 0 I “1 , 3 I 

J 

000 5 


F5=bS+FCI f 1) 

■ 

0006 


K = l + 1 


000 7 


DO K- J~\,32 

• Q ( 

t 

i I « i U ,n 


;w J, •» 1 1 


0009 

10 

F b = F 5 + 2 ^ H 1 1 r J ) / ( 1 + M ) 


0010 


COMK'iUE 


0011 

r5 = F:D + PC 32,32) 


0012 


pF,rUHN 


OOli 


END • * 



I 


FORTRA'J IV VOlB-02 FFM 26-N0V-76 igM-lrSB RAGE 001 

CORr: = OflK, LIIC=l2t>0,200J • vSUBKOl’S , LR : /LI : 1 =D7 1 : SUf-.ATJTS , 


0001 

0002 

0003 

0004 

000!> 20 
0 ‘j 0 o 

000 7 


S U o R 0 u T 1 i\i K 5 A b C F 6 3 

PC32, 32 ) , ALrA(32) ,bLlA(64) ,GAMA(32) 

F6 = 0.0 

DO 20 1=2,64 
P0 = r6*rIrD(i;TMC I ) 

KbIUHN 

ErvD 




FQHTRA^ IV VOlE^-02 * FHl 26- ^UV-76 1 9 : 1 b : 1 7 • PAGE 001 

C0‘^E = 0v 3K, UIC=(200,200] SUBH.OTS , bP : /L 1 : 1 =D 1 1 j i-UbROTS . 


0001 
0002 
000 3 
0004 
•0005 
0006 10 

0007 

0008 


S U B.R 0 U n N E S U M V A lU P 6 , F 7 ) 

COMMON P C 32 , 3 2 ) , ALFA C 32 ) , BETA ( 64 K GAMA ( 32 ) 
F7=^0.0 

DO 30 1=2^64 
H=l-F6 

F7=F7tHf H^BETAC I) 

RETUKN 

END 


FnRTRA(V IV VOlB-02 FRI 26-N0V-76 V9:i5t40 PAGE 001 

COP£=^03<, U 1C =[200,200] SUBKOl 3 , bP t /bl : l=D‘l 1 rbUBHOTvS. 


0001 

0002 

0003 

0004 

000 5 10 

0006 

0007 


SUBROUTINvE SUMENTtfQ) 

COMMON P C 32 , 32 5 f ALFA ( 32 ) ,bETA(64) ,GAMAC32) 
F3=0,0 

DO 10 1=2,64 

PB = F8-BE1 a C I ) tALOG (BETA ( I ) -^ 1 . OF-7 ) 

RETURN 

END* 


FORTRAN IV VOlB-02 

CQRE=U8K, UIC=[200,200] 


FRI 26-NOV-76 19:16:03 PAGE 001 

SUBROTS, LP:/ LI ; 1=D11 : SUBROTS. 


0001 

0002 

0003 

0004 

0005 

0006 
000 7 

0008 10 

0009 20 

0010 
0011 
0012 


SUBROUTINE FNTCF9) 

COMMON PC 32, 32 3 , ALFA (32) , BETA (64) , GAMA (32) 
F9 = 0.0 

DO 20 1=1 ,31 

F9=F9-(P( 1 ,1) )*AL0G(PCI,1)+1 .OE-7) 

K = l.+ l 

DO 10 J=K,32 

F9=F9-2^CP(l,a})^AL0G(PCl, J3+1 .OE-7) 
CONTINUE 

F9 = F9-P( 32, 32) ALOG C P ( 32 , 3 2 ) + 1 . OE-7 ) 

RETURN 

END 


FORTRA'I IV VOlc-02 

CQHK=03K, lJlC=t200,200] 


(•NX OV -7ri IV.'lbSiV 

SUbKO'KS,LN'!/Ll: 1=DT1 


0001 
U 00 2 
00U3 
OOO-l 
0005 10 

OOOo 

0007 

0008 

0009 20 

0010 
0011 


SUhRUllfr.fr: DlFVARfFlO) 

CQmIMU.'v P ( 32, 32) ,ALFA( 32) AC64) ,f':AMA(32) 

UIFAVb, = 0.0 
DO 10 r=v,3i 

D 1 F A V E = D 1 F A V P; + 1 f G A M A ( r + 1 ) 

FlO = u.O 
DO 20 f=0,3l 
H=l-DIFAVe 

F10 = F10 + >i + irGAMA(IH ) 

ReruKiS' 

EbO 


FOBTRAV IV VOlB-02 FHl 26-DiJV-Vb 19:16:53 

COP.E = 08K, U1C=C200,200) 5U6RC1TS , Lt- : /LI : 1 -D’i 1 


0001 
0002 
000 5 
■ 0 0 0 4 
' 0005 
,0006 10 
( n 007 
. OO'^ei 
\ 


•SUOROiJTtl.’E DiFEIvTCFl 1 ) 

PI j2,32)fALFA(32) ,BtlA(64) ,GAF.ACJ2) 

P 1 1 - 0.0 
DO 10 1=0,31 
n/sGA^AC 1 + 1 ) 

F11=F1 1 + ALOGCW-H .0E>7) 

Pi. ruPM 

t.f,U 


FORTRAN IV VOlb-02 Ff 1 26-NDV-75 19:17:18 

CORE=08K, U1C=C200,200] ■ SUBHOTS , LP : /LI : 1 =DT 1 


0001 


SUBROUTINE C0P.r.E(F9,FJ 2) 

0002 


C0..1H0N P( 32,32) ,ALFAC32) 

0003 


HZ=0.0 

0004 


Hxyj=o,o 

0005 


DO 20 1=1 , 32 

0 00.6 

* 

H = ALFA(1)t-ALFA(I) 

0007 


HZ = .-!Z-(ALFAC X) r=ALOG(ALF. 

000^ 


* = 6l,C:i:(H+ ! .uE-7) 

0009 , 


hxri=riXYl-A*P(l,I). 

0010 

• 

P = l + 1 

0011 


IN CM.GT. 32 )GOTO 30 

0013 


on 20 J=M ,32 

0014 


F=CALFACI) )tCALFA(J) ) 

0015 


C = ALrjG( h+l .Ot-7) 

0 016 

2 0 

H/Yl=rt;\Vl-2 + C + FCI,J) 

0017 

30 

n /sCN 9-HXYl )/HZ 

UOId 


HElUHU 

0019 


END 


• 7 ) 


PAGE 001 
:SLhRDT.S. 


PAGE 001 
: SuBHO’IS . 


PAGE 001 
SU^POTS. 


B.13 


]?AGP to 

POOE 


FOKTRAV IV V013-02 Khl 26-|vOV-7& 19:17s44 

CORS = «)iiK, U1C=1200,200J SUBUGl’S , LiF t /LI ; 1 =D‘i 1 


0001 SUHROUTIkE FQXT(VW) 

0002 D1’-k.;.S.I(j!'! 12) 

0u03 CO R ( 32 , 32 ) > ALFA ( 32 ) , F Kl A( 64 ) f GA.'iA C 32 ) 


0004 

CALL 

AUX 

000 5 

CALb 

AUXl 

OOOo 

CAM. 

AUXO 

0 0 0 7 

GAbL 

ASMCV.vU ) ) 

000^ 

CAL’u 

C0M(V'.C2)) 

0000 

CALL 

CUHCU,Vrt(3)) 

OulO 

CAbL 

VAR( U, V.V( 4) ) 

0011 

CALL 

XlD^i('/.v(5)) 

0012 

CALL 

SL -1AVE(V.;(6) ) 

0013 

CALL 

RU-'.V AR(VW(6) , V.v(7) ) 

0014 

C A.L L 

SilMENUVrt(8) ) 

0015 

CALL 

thfC V‘Vl9) ) 

0016' 

CALL 

Lllf VAR(Vw(10)) 

0017 

CALL 

DlFSIU-l VV. (1 1 ) ) 

0 0-16 

CALL 

CUKMt(V^.(9) , Vi‘K12)) 

0019 

HetURL 

0 0 > i ) 

• LNLi 



RAGE 001 
bUBRDIS. 


KOfiTi^AN IV VOIB-'O? 

eUKK=03K, IUC=(200*200] 


hhl' 26-NOV-7>j 20:26:10 PAGE 001 

MGKAH.I.P: /l,J ; 1=DT1 tHGRAM. 


OCOl 

0002 

OOOJ 

0004 

0005 

0006 
0007 
0003 

0009 

0010 
0011 
0012 ' 
0013 
nOl 4 

dois 

0017 
V.OX’6 
no 19 
0 021 
0022 
00 2 3 
002 4 
/ 0 A 3 
0026 
0028 

0029 

0030 
00 31 
0032 


0034 

0035 
00 36 
0 0 37 

0038 

0039 

0040 

0041 
0042. 

0043 

0044 

0045 
0 0 46 

0047 

0048 

0049 

0050 
0 05 2 
0053 
00^4 
0056 


C 

C PKOGRAMA PAIU V I SUAL. EZACAC) HA CAE^ACIDAOL DK SLPARACAO tNTRfe: 
C CUf\$SES USA^^5Da CARACrphlSnCAS Tb\iuRAli> LVOU KSPEC7hA15 


IfH'EGER SCC2/2). 

LQGJCAL^l vU7 4'0 

DlSEN.SlOtv X(20,200) , bX ( 200 ) , COU/‘. I ( 50 , 2 ) / C 2 J 
DIMb^SXO^i PEAl (44) , YrAX(2) , ICLA(2) 

OATA lY/440,0/ 

CAbb OOTPUT(27,12) 

Wf?lTeC6rn 

I FORf-iATC l3Xf ‘ 44^: SEb6*CA0 DE CARACYbB USTICAS 

5 aRITE(6,10) 

10 fORKATt * saUAKlAS AMOSTRAS EX IS TEM? C 1 -200 ) > » ) 

CAbb 0UrPUT(7) 

READ ( 6,20 }iv 
20 raRViAT(74Al) 

CAbb FROMT(^,74) 

IF(vyO) »E‘0.*‘i30) GO TO 2000 
•b=0 " 

CALL IN-TFr(L,«,74,f;L’.v) 

Ir fi»UM,LT.i .nP.L'lJW.GT.JOO) GO TO 6 
25 wKI i'E( 6, 30) 

30 t'OBNATC ' SFA1XA(MAX = 20) DE CARAC I'ERISTICAS : # IMCIAL E F1NAL>') 
CALL 0UTPUT(7) 

HEAUfb.2’Jlisi 

C A L i j ? ?•>.(•. I Lfirl i j 

IFCivC n . EQ. "1 30) GO TO 2000 

L=0 

CALL 1NTFF(L,W,74,MIN) 

GALL INlFF(L,'w,74,MAX) 
lFAX=f'«AX-Mllvi-l 

1F(MAX.GT.44.0R.M1N.LT. I .OP.MAX.LE.I1IN,OH.IFAX.GT.20)GOTO 25 ‘ 
C LEIUHA’ DOS DADOS OUf INTEHESSAM .70 AROUIVO TAB. DAT 
C DA DECJAPE^d CRIADO PEXO PFOGHAMA 'COLETA' 

CALL FDBSLT ( 3 , ' READONLY ' ) 

CALL ASS1GNC3, ' D1 0 ; TAB . DAT ' ) 

DEFINE FILE 3 ( 200 , 1 00 , U , NREC I 
MREC=1 

DO 39 1=1, NUM 

HEAD C 3 ' WPEC) ( f ( FEATOJ ) ) ,-J= 1 , 4 4 ) , SC , LX (1 ) ) 

DO 31 L=MIN,MAX 
'•■ = L“HIN + 1 

31 XC''., l)=fiEAT(L) 

3^ COmIIWUE 

CALL CLOSECi) 


40 

50 


wPIIE(A,50) 

FOB.SATC ' SDDE CARACTERISTIC A NA FAIXA 
CALL 0UrPUT(7) 

HEAD(d, 20)W 
CALL FRQNT(w,74) 
iF(NCl).EO."130) GO TO- 25 
L=0 


ESCCLHIDA GLEP DSAR?>') 


CALL I.'fTFF(L,H,7 4, IFFAT) 

IFUFEAI .LT.«IN.OH.Ire AT.GT,MAX)GO TO 40 
uFE AT = I FE.AT“ M I Ni l 


FC'»rHAN IV VOlB-02 

COMKsOBK., 0IC = tV.OO,2O0] 


ft'I ?t•'^.0V“7o ?o;26;i0 • PAGE 'it/ 

, HGKA;>5,1 P.V/Ll : J =D'I l lUGrA*'. 


1 


OQ‘5*; 55 

0058 50 • 

0059 
OOoO 
■0061 
00b2 
0064 


wKU't'H.eo) 

Ki'JO’>'AT( ' St'SCQLHA UM E AF LE CbASSflS ( 1 ~b ) > ' ) 
CALL 0UTPUU7) 

Rt.ADCb,20)W 
CALL FHDNr(.\, 74) 

It'(W(l).fc.U.")30)Gn fO 40 
L=0 


0065 
OOob 
0067 • 

0069 

0070 
007 1 

0073 

0074 

0075 


00 68 1 = 1,2 

CALL lNTEF(L,'.-,,74, ICLA(I) ) 

JFUCLAll ) .G1 .O.GR. ICLAU ) .LT.l ) GO *10 55 
68 CONTINUE 

C DIvSCKE:T1ZAC;AD las DISTRIEjUICOLS E ObTE.'.CAO DOS HISIOGI'-A.MAS 
DO 100.1 = 1,IJUF 

lE'ILxm.OL.lCLACn.AND.LXm.NE.ICLAIZllLO TO 100 
XViAXsXIJKEA'I , 1 ) 

H?UN=XhAX 
llviC = l‘H 


0076 

00 7 7 10-0 

C078 140 

00 7 9 


GO TO 140 
CONTINUE ^ 

DO 200 I = IMC,KUM 

U' ( LX 1 1 1 . Nt . J Cl-A ( 1 ) . AND . LX n 1 . ME . 1C La ( 2 )) GO 'lU 20 0 


00 81 
Q0«2 

0083 20'0 

0084 
0 0 


X^^AX=A.MAXH;a JKEAT, I) ,XPAX) 
= { X(iiFK AT, n ,XM1N) 

COM InUE 

DEL'^A=CX^'AX-X.MIN)/48 
D U 23 0 J 1 f 2 


uu bb 
OOd? 
0080 
0.0 8 9 

0091 

0092 
009 3 
0094 

0096 

0097 

0098 

0099 

0100 
0101 


l)U /5U 1 = 1,50 1 
250 CDUMd, J)=0.0 
DO’ 500 1 = 1 ,NUM . 

IF(LX(n.NE.lCLA.(l))GO TO 300 
J = JFIX((X( JFF.AT, I)-XLlfn/LELTA)-H 
CiJUNfCJ, 1 )=COUNTC J , 1 ) + l . • 

GU TO 500 •• ' 

300 lF(LX(r).N£.lCLA(2))GD TO 500 

J = IF IX ( (XUFEAT , n - XLIN) /DELTA ) + l 
CDUNT(J,2)=C0UNTCJ,2) + 1 . 

500 CONTINUE 

C CONSTRUCAO DOS HISTQGRAMAS NO TERMINAL GRAFICO 
'JO 600 1 = 1 ,2 
VMAX(1)=0.0 
00 590 J = 1 ,50 


0102 590 

0103 500 
010 4 

0105 

0106 
010 7 

0 i 0 8 

0109 

0110 

0111 550 

0112 

01 13' 

0114 

0 1 1 5 700 


VNAX ( I ) =AMAX1 C YMAX C 1 ) ,CDUNTC J , I ) ) 

YMAXm = l .inMAXCI ) 

CALL I.NITT 
DO 700 1=1,2 

CALL UaInDOCO. ,51 . , 0 . 0 , AX ( I ) ) 

CALL S'.-;1NUO(10, 1010, lYCl) ,340) 

DU 550 0=3,50 | 

ALC=FLOAT(J) 

CALL MOVEA( AbC ,0.0) 

CALL DRA«A( A6C/CUU.‘.'1CJ,1)) 

CALL. «OVEA(u.O,YMAX(I)) 

CALL DrAWA(0.0,0.0) 

CALL DRA.4A(51 . ,0.0) 
continue 


0115 • CALL MUVABSt O , 430 ) 


FOKTRA-J IV 
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20:26; 10 

l’Ac:c OOi) 



0 ^ K y 

ulC=t?.UO,200] 

h i, P A M , ij P : / L J : ' 

1 ^ 01 1 J Ai** • 


on? 

* 

CALL A'iMQLE 



i.- 

0 i 1 


rtKll bCo , bOCO U LAl ,1CLA 




0 n ^ 

’^00 

rUff'ii'! Cl X, 'CAf-ACrtKlSTlCA liSAD.A: ' , 

13, 5X, 'CLAsSSXS 

Ml. 




C ' !■: ',11) 




01 ?Q 


hKAi'lo , 20) >v 




012.1 


CAl b Uu*'lPUrC27, 12) 




0122 


, GU 10 ‘'S 



■■■ * 

0123 

2000 

STOP 




0 1 2 4 


E I'v D 





- 31 - 

Cada janela foi equal izada em 32 nTveis e colocada no C£ 
nal 5 do si sterna I-lOO, de onde eram extraTdas as respectivas matrizes. 
Como foi visto no CapTtulo III, cada tipo de textura fornece urn conjunto 
de 3 caracterTsticas, que sao a media, o desvio padrao e a faixadevari^ 
?ao dos valores dessa medida, para as quatro matrizes de coocorrencia. 
Tem-se, entao, urn total de 36 medidas de textura, as quais serao denom_i_ 
nadas pelo tipo de medida seguido da letra M, D ou F, caso se trate, re£ 
pectivamente, da media, do desvio padrao ou da faixa de variagao. 

Alem dess as 36 medidas de textura, foram, tambem, compute 
das 8 medidas espectrais: media e variancia dos nTveis de cinza dos poji 
tos da janela, para cada urn dos 4 canais, perfazendo urn total de44 cara£ 
terTsticas para cada janela. 0 tempo total de processamento, de cada j£ 
nela, foi por volta de 40 segundos. 

Na fase de cl assificagao, o conjunto das 182 janelas foi 
dividido em duas partes iguais: uma para treinamento do classificador e 
outra para teste. Portanto, o numero maximo de acertos nas label as, que 
serao mostradas mais adiante e 91. 

4.3 - RESULTADOS 


Inicialmente foram realizadas classificagoes , utilizando 
cada uma das 44 caracterTsticas i ndi vi dual men te. A Tabela IV.l mostra os 
resultados obtidos, com o uso do WTn-max e do vizinho mais proximo (VMP). 
Nesta Tabela, as caracterTsticas espectrais sao designadas pelonumeroda 
banda correspondente (4-7), junto com a letra M, no caso do ser media,ou 
V, se for variancia. 

t difTcil fazer uma analise baseada apenas nos resultados 
de classificagoes, usando caracterTsticas individuais. E, porem, imprati_ 
cavel se fazer todas as combinagoes possTveis de caracterTsticas, para 
analisar, em bases soli das, a importancia de cada uma del as. Num proce£ 
so natural e empTrico de escolha, pode-se selecionar aquelas que deramm£ 
Ihores resultados individualmente, desprezando-se as demais . Ou entao, pr£ 
cura-se uma segunda caracterTsticas que, junto com a melhor delas,indivi_ 


FOPTRAIJ IV VOtH -02 

CURb'sOflK. UIC=( 200 , 200 J 1 


FRl 26 -fiO ’*-76 20 ? 26 :l 0 PAGt' 001 

CAIEG,I*P:/Li; ls:DTl tCATKG. 


0001 
*0002 
000 3 
000 4 

0005 

0006 

0007 

000 3 
000-9 
0010 

001 I 

0013 

0014 

0015 

0016 
001b 
UOiy 

0020 

0021 

0022 

.0023 

0025 

0026 
0027 
002d 

0030 

0031 

0032 

0033 

0034 

0035 

0036 
003 S 
0039 
0 040 

0042 
004 3 

0044 

0045 

0046 

0047 

0048 

004 9 

005 q 
0051 


w 

C PROGP.AMA PARA CbASSIFIC ACAO A‘'.OS|l<AS 

C POi>EM SER USADOS 2 vf.’TOCHuS : 0 M 1 .\-VAX E U CO VI'^INHO MAIii PhOXlPO 

c cfir.'jucro de. thli'sci eh cousti iuido oa.s apostpas ikpares 

C COIvJUN'TO UE lESTE EH CONST I TU 1 00 DAS AMOSIRAS PARES 


V 

LOGICAL*! WC 74 ) 

DIMENSrOH 1 XC 2 . 2 ) , JCOMPC 20 ) , liVtC ( 5 ) , LX ( 2 U 0 ) .A( 20 , 5 ) 
DlMENSin.N 8 ( 20 , 5 ) , )‘( 20 ) , PEATC 4 4 ) , X C 20 , 200 ) , i’iUT C 3 ) • 

CALL' OUTPUTi 27 , 12 ) 

1 WRITE( 6 , 2 ) 

2 fORMATCl 3 X,'t** CLASSIFICADOB ***',//,' TOTAL DE', 

C AMOSrHAS(l- 200 ) E DIKEHSAU 1)0 ESPACO ( 1 -20 s> ' ) 

CALL Ol)rPUTC 7 ) 
h'EAD(6,5)W 
5 FORFATC 74 A 1 ) 


1 3 
I S 


20 

21 

22 


CALL t Runt (IV, 74 ) 


nK.v(l).E(J."130)GO TO 2000 
L=0 

CALL i^ri Tf (L-, w, 74 , ^'uPP 3 
CALL 1 NJ'FF(L,W, 74 , JOIM) 

IF( rJlJpiB. L E . 1 . OR . NUMB . GT . 200 . OR . JO J '-5. 
/■R1TE(6, 15) JDIM 

Fnpr/iTC 1 t , 'E*-TPE CG ■ AS CA.'iAL 


(:T. 20 ,OR.JD 1 M.LT .1 )G 0 TO 1 


^ 4. U A . 


wiJD Vhl t 


V r / / " 'o ^ ' J 


CALL 0 UTPUK 7 ) 

READ(( 5 , 5 )W 
CALL FRONT! W, 74 ) 
iF(W(l).£Q." 130 )GO TO . 2000 
L =0 


DO 20 1 = 1 ,J 0 IM 

CALL' IN J'FF ( L , W , 7 4 , JCUMP ( 1 ) ) 

IF(JCOMP( 1 ).EQ.O.OP..JCO(.'P(I).GT. 44 )GO TO 13 

CONTI HUE 

i'.BirE(0,22) 

FORFAl ( ' S JUAN'TAS CLASSES EXl STEP ? ( 2 -b ) > ' ) 
CALL 0 UrPUT( 7 ) 

READ ( 6 , 5 ) w 
CALL FROlvT(w, 74 ) 

IF(iv(l).£Q." 130 )G 0 TO 2 ()C 0 


L =0 


CALL I\TF.F(L,^;, 74 ,KCLA) 

• 1 F(KCLA.LL '.1 .DR.IvCLA.GT. 5 )G() TO 21 
C LEI TUBA OF. DADOS £ FONTAGEM DA HATPIZ DE AMOSTRAS 
CALL Ff)BSEI( 3 , 'READONLY') 

CALL ASS 1 GN( 3 , 'DKOrTAB.DA) ' ) 

DEFlIvE FILE 3 ( 200 , 1 OQ , U , KREC ) 


tv«EC=l 


DO 100 1 = 1 ,NUMB 

READ( 3 ' NBEC)( ( (F EATC J) ) ; J= 1 , 44 ) ,IX, LX(I) ) 

DO 30 J= 1 ,JDIM 

K=JCOMP(J) 

20 X(J, 1 )=FEAT(K) 

100 continue 


P0RC«A1 IV 

VOlli-02 FM 2b-NOV-76 20:2o:l0 PACK 002 

CORES 

' J 8 .< ♦ 

UIC=1200,200] . CAfFG,LP:/L,l: l=DI’l :CArF.G. 

0052 

■ 

CALL CLOSiLU) 


C NfiRMAblZACAU DOS FKATURKS CKN'JPCI Ot CAOA COMPCmKME, 

0u53 


UC 210U J=1,J01M 

0054 


X.MAX = XU, 1) 

005‘i 


XX 1 NSXVAX 

0056 


Ul) ,iOSi) J = 2,NUF13 

0057 


xmax=:a..;axi cx*'AX,X( I, j) ) 

0055 

2050 

X’-U A=A-.- in ( X^II S , X ( I , J ) ) 

0059 


DO 2UH0 Jsl.NUFH 

0 0 6 u 

V060 

X(I , J) = ( x( I , J)~XMIN)/(XMAX-XM1N) 

; 0061 

2100 

CfJNTl rjllfc' 

0062 

104 

tvAlTLlb, 105) 

0063 

105 

FOPMA'H ' SOUER USAR 0 (M)INMAX Oil 0 (V)lZlNhu MAIS Pf UX 1 vt.i ?> • ) 

0064 


GALL UUIPUTI?) 

0065 


Hb-AU(5,b)rt 

0066 


CALL FRC'‘a(W,74) 

006 7 


lF{.v(l) .En"l 30)G0 TO 2000 

0059 


JKw(1).i-0."12is)GO TO 530 

0071 


IF(.,(n .It. " 1 15)GO TO 104 


C CLAGSXriCACAD COM 0 METODO M3NMAX 

ir% 

0073 


UO 110 1=1, KC LA 

007 4 

110 

I NIC (I) = 2 

0075 


DO 200 1=1, NUMB, 2 

' > U 7 0 


j=Lxm 

» - 
V V < < 


I. lu cOv 

007 3 


DO 150 ri=l,JDIM 

008D 


A ( W , J ) =X ( M , I ) 

OOHl 

150 

BIM, J)=XC<1, 1) 

00tJ2 


IiaC(J) = l 

0003 


00 17U N=1,KCLA' 

00 B 4 


IF (Iwicm .N£. 1 )G!1 ro 200 

0085 

170 

CONTINUE 

0087 


GO TO- 250 

0088 

2 00 

. eONTlfiilF 

0«B9 


tf<nE(6,210) 

0090 

210 

FOPMA'l ( lx, ' ENGAND-NAO EXlSTEN-l TANTAS CLASSE51' ) 

|0091 


GO 10 21 


DO PAHALELEuFIPl-.DO CUFREGPOfvDEFTL A CADA 


I ' C CLASSL USA^uD U COpJUino DE TFElND( AM03THAS 1FPAHE3) 
|U0^2 250 DO 220 <^=1 ,KCLA 


I 000 3 220 

|0004 
'0095 
’0096 
0097 
009 5 

0099 230 

0100 400 

0101 
0102 

0103 

0104 

0105 430 

0106 


1 viC(N)=0 
DO 4 00 I = 

N=DXCi) 

j01Clfv)=l^iIC(^O‘+l 

L;0 230 v.= i,JDiM 

A( v,rO=:A* AXl CA(M,N) ) 

h ( , w ) =A V, 1 N 1 ( b ( M , N ) , X ( AM ) ) 

ZOuriHLl 

DO 430 J = 1 ,KCt;A 

DO 430 , JOIM 

T E r A - ( A ( J ) - B t M , J ) ) / C I N 1 C ( J ) - 1 ) 
J)rA(.V;, j)+tETA 

b / i* 1 \ ^ , 4 / , I » % _ »Tt 0 /r » • 

V . J , SJ ; - U V • 0 \j j I tL I H 

DO 440 hOD=l,3 • 


|i 


'A: 




! 

* ^ p 

CORK = 

08K,. 

& ' 

1 i 

J 

0107 

440 

j 

•( 



C TF- 

i 


0108 


f . 


0109 


?! P 

t' ■ 

i' 

0110 


f 


0 m 


i 

0112' 


1 

Oil 3 


1 ; r 

0114 


1 ■■ 

01 1 5 

*470 

1 

0116 


i 

-0117 


i 


C CO 

\ ' ' 


C A 


01 18 


! - 

0 120 


■i 

: 

0121 


i 

0122 

480 

I 

0123 


T 

1 ■ 

0125 



01 2D 

490 


0127 

son 


0128 ’ 



0120 



C \ 3 

\i ’ 0 

1 

0 1 j i 

» 

II 

.0132 



0133 


V 

0134 

520 

i 

0135 



0136 

521 


0137 



0138 

522 

. [ : 

0139 



0140 

525 


VO 1 1^-02 


FHI 2b-f CV-76 


20;?6;i0 ♦ PAGE 003 

C ATEG , LP : / L 1 : 1 r-i)T I ; C ATEG . 


0 1 ‘H 
0142 


'E in: CAOA UMA 
1AM = 0 


D/S AMOtSTFAS PAHF-S 


IEHH=U 
IQA=0 
XKK = 0 
DO boo 


1=2, NUMB, 2 
lAMslAMI 1 
DO ^ ;0 J=1 , JDIM 
V(J)=X(J,I) 

CAUL TESTA ( J D 1 M , KC LA , A , L , V , IC AT , MOO ) 
iTor(:.ujCO = iroT(MGD) + i 
MPAKA A CATEGOEIA ICAT ATRlbLIOA P£LQ 
CArPJGORlA VERDADEIRA OUE: ES*1 A CfuM IDA 
1KCJCAI.NE-LX(1))G0 TO 480 
lUK=IOKtl 
GO ro SOO 
IERR=:ILHRf 1 
IFCMOU.LQ.i) lEK=lLR-f 1 
WRl ir;Cc.490)LXCI) , ICAT 
FOk^aU IX A *0^3 ISA Oh CAIhGUt-lA 
Cr ‘TEGOiVlA; M3) 


CLASSlFiCADOR 
NO VETOH LX 


COM 


12,' CLAStslt iCAU« COXO C4' 


CRN l‘l 

PEPCbN=l j^^lOO.O/IAM 
wRlTE:(6,blU)PEPCEN 


r •- ,* c i T '''* *- r * 


i l *.6 n- - ij t* - i A 

Call assignc 5, » lp: ' ) 

V. U r 1 1 ( b , b 2 0 )'( tv TK ST£ , 1 AM ) 

F0RMAT(///,10X, HLSULl'AUO DA CLASiilt ICACAD-MIK'MAX »♦•»=' 
C3X , ' ASOSTRAS PARA rPtl NO : ' , 1 3 , 5 X , ’ AMQSTR AS 1 bSTAOAS : ' , 13 ) 
V.'PlTb(b, 52U ( J01M,IERR, ( (aCOMPCJ) ) , J= 1 , JOlK ) ) 

FURWA'l (IX, 'UliSt-NSAO L'O t SPACO : ' , I 3 , 1 OX , ' J L TAO DL', 

C fRPDS: ' ,13,y, ' CARAC1LP1S7ICAS: ' ,12I3,/,17X,8L3) 
vStUTtCb, o22)PERCEN 

KORHAKbX, 'Pr-RCtNTAGEM DE CLASSif 1C ACAO C'jRRLT A : ' , f b . 2 , ’ % ' 
7.RllE(b,b25nTOT,lER 

FORMATU'*, IX, , 'CAIRAM NOS blMlTES DE >^PEivAS UMA CLA5SL«,/, 
CIA, IX, 'CMRA’-‘ NT'S LlKllES HE KAIS DE l!J‘A CLASSK ' , / , 1 4 , 
ClX,'NAu CAIPAf-' DEKlRtJ DCS L-l'-'inES OE ALGUl-'.A CLASSE', 


C/, J4,1X, 'DESIAS 
CAL'!, CLUSECS) 

GO TO 104 


ULl’ll-AS RESuLIARAM cM LRRD ' ) 


C CLASSIFICACAO BA5EADA -NA CATEGORIA DO VIZINHC KAJ.S PBOXIKO 


0143 530 

NTLSTE^O 

0144 

Id^=o 

0145 

ie:rr=o 

.0146 

DU 9UU 1=2 r NUMB, 2 

0147 

N1LSTF=» I'ESTE+l 

014 8 

■ DO 790 J=l , NUMB, 2 

0149' 

PRBM=0.0 

0150 

DO 7bO K=1,J0JM 

Ol5 1 

j = X 1 ^ , U - A ( r # J ) 

0152 

l = Ati6C 1 ) 




fix... 








1' tjR'iM 

\*i IV 

C0KE = 

03K, 

01 S3 
0154 

0 1 S o 

0157 

780 

0159 

0160 

0161 
‘0163 
01 bi 

;90 

C CO 
C A 

0165 

0166 

880 

0167 
0163 
0169 
017 « 
0 1 7 1 

900 

017 2 


0 1 7 3 

910 

017 4 
01 7 5 
0176 

1 7 -7 


‘ J 1 / 7 

0179 

•) u u 


Ill 


Vu1l<-02 

:-l200,200\ 


Vhl Vj6-l-0V-/b .^0:2n:10 PACK «wr04 

CATSt.^'r LPl /Ll: l=DTI :CAft:G. 


PKLu CLAS^It ICAilOH 
liUA IvO VETUH LX 


PRBM = PP^M+‘1 
If ( J.EO. 1 )5tLl=:PPbM 
Sn.LI-A'!PcJ CPAf^f^-,SELn 
1 f ( SLL I . to . PHBM ) 1 N'OtX =v3 
CUKriNUE 
1 ChT = I.XC IM)tX) 

ypApA A CATEGuA'IA ICAf A'UUbUIDA 
CATtGORIA VFPOADtIRa uuf r..STA COL 
irClCAT.tvF,LXC 1) )GO 10 ft80 
iUK=iOK+l 
CO TO 900 
jEKR=lFRR+l 

AKITF-Cb, 190) (LXCI) , ICAT) 

CO’NTINUE 
I AM=NUMb-MrtSTE 
?tKCEL = i:iK4=l00.0/IA.M 
/.•h'lTECo, biOlPtPCEN 
CALL ASSIGLC:), ' LP; ’ ) 
i 1 1 1 :3 , bi; 0 ) C H i tSTE r 1 AK) 

^CRMAT(///,2a, ’ RE SIJLTAUO da CLASSIFlCACAu-VIZIi 


COM 


C/ 


kBOa 1 ^ ^ 


PAPA iHhi^D: 


, iX / ‘ Al' Ui? 1 KA5 
/ * j itb • , ii ) 

v^hiibCS^b/l) 1 JDIM, lEHP, { (OCGMPCOn m1=1/OOTM) ) 

;vKnE(b,52?)PEPCEL 

CALL CLOSE(b) 

0^ it iv-1 
o 1 u r' 

\^H0 


13. 


. H 0 K’ A I .‘P 
. » A:a:ks ‘ 


1 

T 

I 

T 




FOPPRAM IV 

CORK = OtiK, UiCti L'.!U0,200] 


0001 

0002 
QQua 
OOO'i 
0005 

OOOt? 

0007 
0 0 0 0 

00 I 0 

001 1 
U 0 I 2 

-0013 
i 0015 
I 0017 
001-5 
0020 
00 21 
0022 

0023 

0024 

0025 
002 6 
002 H 
0029 
00 30 
0 0 ^ 1 
u\j J /i 

. 0033 
0034 
0036 
0036 
003^ 

0040 

0041 
00 42 

( 0043 
i 0044 
: 0045 
0 0 4 h 
0043 
0 049 

0050 

0051 
0052- 

0053 

0054 

ooss 

0057 

0058 
0 05 9 
0061 
0062 
0063* 
0064 
0 0 tj 5 


H<1 2r-Nnv-/6 20 : 2o:l0 . PA(Jr: OOJ 

CAThG^LPs/LI t l=iDTl ;CArbG, 


SUBRUUTircK TF81AC.JDn.;,KClA/ A,1W Y, JCAT.MOU) 
OXMfc:^vSUn A(20 ,5) , -^(20,5) , t t 20) ,FR5(5 ) ,L1I C5) 
•'iU^O 

DO 15 l=:l ,bCbA 

0I.4C n=o 

00 5 J=1 , JDIF 

IFCYCJ).f/T*fHJrI),OR.YCJ).GT.ACJ,n)GO TU lb 
5 ca.vnwub: 

LIMC I}=:1 
NO=NOf I 

15 G0^^ ] IK'JiL 

lK^vn,EO,03G0 10 500 
lF(*NO.Gh*i )G0 TO 300 
DO 30 1=:1,^CL.A 
IF(LIMCI) . b'Q. 1 )G0 fO 40 
30 CONTINUF 

40 ICAT=1 

HOD-1 
KbiuRN 

300 PRBH^O.O 

i;ij 35u J = I , KCi.H 

1 F 1 L 1 ''M .[ } . b ,\ )irn '! g 350 
PKbCl)=^l .0 

00 330 jri , JDIM 
3 30 PRH(l)=pi^JKI)/CA(J, 

> vr:. 1:! \ c r--. >*. f i } ^ -v-Rv. } 

J.6 b C u»> i a u u b , 

UO 360 1 = 1 ,KCl-A 

lF(Li;-'(i) .i.tr. 1 )GO TO 360 

.1FCf-His(U .c;O.PK^J^l)Gu 10 380 
360 CCiM’lNUfc. 

380 1CAT=.1 

MOD=2 
KETURN 

500 PRB(-i=0.0 

t)0 600 J = l,rtCLA 

PRbC0)=0.0 

DO 550 1 = 1 , JDIF. 

IFlYfllcCM .c'(I,-J).QH.V(I).LT.A(l,J))GD TO 560 
rETA = i«, 1)-B(] ,J) 

TEl A = AciSCJ ETA) 

VE1"A=U1)-A(1,J) 

VETA=Af'S( VETA) 

TETA-ANai'il CVt:iA,TETA) 

PHnCJ)- ^B(J)4lETA^rETA 
550 CONf' 

1F( , ,0. n PR5M=PPf3>..- : 

500 ?RBM=A n‘;t CFP£f-'.,PF6.J)) 

UO 710 J=1 ,KCLA 
lF(Pnbi-;.Ei).PRt>(J) )GO TO 720 
710 COixTlNUE 

720 ICAT=J 

(•100=3' 

RETURM 

r-* I r> 




I 

I 






C IV U 


